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ABSTRACT 


The minority carrier diffusion length is a critical 








parameter in th development of next generation 
Heterostructure Bipolar Transistors (HBT) and highly 


efficient solar cells. A novel technique has been developed 





utilizing direct imaging of electron/hole recombination via 
an optical microscope and a high sensitivity charge coupled 


device coupled to a scanning electron microscope to capture 





spatial information about the transport behavior (diffusion 


lengths/drift lengths) in luminescent solid state 








materials. In this work, a numerical model was developed to 





do a multi-parameter least squares analysis of transport 
images. Results were applied to the study of transport in 
materials at the forefront of device technology that are 


affected by quantum scattering effects, where few reliable 








experimental measurements exist. The technique allows for 








easy localization of th measurement site, broad 
application to a range of materials and potential 


industrial automation to aid the development of high speed 





electronics for terahertz devices. 
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I. INTRODUCTION 


A. HISTORY 


With th mergence of the transistor in 1947 came a 





revolution in military affairs (RMA) that has been evolving 





over the subsequent 59 years. Today, the battlefield and 





our daily lives are littered with electronic devices that 
do everything from helping us to see in the dark to 


enabling global communication links. The transistors that 





act at the foundation of these capabilities are able to 
perform faster and faster every “Moore cycle”. As the 
demand for faster processing in smaller electronics 
packages has grown, electronics makers have turned to a 
class of transistor called the Heterojunction Bipolar 


Transistor (HBT). 





HBTs are transistors in which at least one of the two 





transistor interfaces is formed of two distinct materials 


[1]. The primary advantage of HBTs is their greater 








mitter fficiency, defined as the ratio of current 


injected into the emitter from an external source to the 





leakage current from the base to the emitter under active 








operation. This advantage results directly from the 
valence band discontinuity at th mitter-bas 
heterojunction. The larger barrier to minority carrier 





injection from base to emitter allows for a substantial 








increase in the permissible doping level of the base layer 


of the HBT, which reduces sheet resistance and allows for 








thinner base layers without concern of mitter-collector 





leakage in the cutoff mode of operation. These advantages 
result in a faster base transit time, defined as the time 


required for the emitter injected carriers (minority 





carriers in the base) to diffuse across the base layer to 
the collector, and a faster switching speed for the HBT. A 
schematic and energy band diagram is shown in Figure 1 
highlighting the valence and conduction band 


discontinuities. 


EERE 
n* — GaAs 






Emitter 





n— AlGaAs 








Base 
Sub — collector Collector 
Semiinsulating GaAs substrate 
(a) (b) 
Figure 1. (a) Schematic cross section of an HBT structure. 


(b) Energy band diagram of a HBT operated under active 
mode. (From Ref.[1]) 


The use of different materials to provide these 
junctions, while adding complexity to the design, adds 
significant power amplification benefits and switching 
speed advantages [.3e]ixs Figure 2 shows a history of 


progression of highest cutoff frequency (F,) for various 


transistors over time. The cutoff frequency is defined as 





the frequency where current gain of one is achieved. It is 


noted that the maximum working frequency of these devices, 


F 


max / 


is the frequency where power gain becomes unity and is 


below F, [4]. 


BOO 





[ —a— EY HEMT 
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Pigure-2.. (Trend) of best. achieved Ff. tor various 
Eransistors.. (From Ref [41) 








In order to achieve thes frequencies of operation, 
manufacturers increasingly rely upon thinner, more heavily 


doped materials to propagate charge quickly and efficiently 





across the base to activate the transistor [1],[2],[5]. 


These increases in doping level and the shrinking of 








relative dimensions, in particular of the bas layer, have 
coupled to bring about new and interesting regimes that 
operate on the edge of known macroscopically determined 


semiconductor transport properties. In order to 











ffectively design and build the most efficient devices in 


these new highly doped low dimensional regimes, new 





techniques that can extract and model material properties 


on the sub-micrometer scale will be necessary [9]. 


B. MILITARY RELEVANCE 


High speed electronics are important for a variety of 





military applications. The Terahertz frequency band, 
defined as the frequency range between 300 GHZ and three 


terahertz, is being explored for applications in medical 





diagnostic imaging, security imaging, and high bandwidth 





communications, just to name a few. A quote from a recent 











Defense Advanced Research Projects Agency (DARPA) Broad 
Agency Announcement (BAA) highlights what might be 
considered the most critical of these applications for our 


information centric battlespace. 


The continuing need by DoD for ultra-high 
bandwidth communications and sensing will require 
electronics that operate at THz frequencies. 
Given that advanced microwave satellite 
communication systems already operate near 60 
GHz, the instantaneous bandwidth required to 
fully monitor the battlespace will certainly 
exceed 300 GHz early in the next century. All 
communication bottlenecks must be removed so that 
surveillance systems can relay their wideband 
measurements to other locations for real-time 
analysis [10]. 

















Conventional electronic sources and receivers are 
limited by resistances, capacitances, and transit times, 
resulting in a significant attenuation of high frequency 
power. High power amplifiers based on HBTs are beginning 
to approach the realm of terahertz oscillations and may 
provide a simple semiconductor based solution for an 


integrated coherent terahertz source and detector. A 





greater understanding of the physics of electrical carrier 





transport in these highly doped, low dimensional structures 





will aid th development of thes devices and provide 








manufacturers with more accurate models and predictive 





design tools. 


Cc. THESIS OVERVIEW 





In this thesis, (sponsored by National Science 





Foundation DMR 0203397) an application of a new technique 
for imaging charge transport is discussed [6-8]. A study 
of a series of low dimensional, heavily doped AlGaAs/GaAs 


heterostructures was conducted with an emphasis on the 





determination of the diffusion length of the minority 








carriers as a function of impurity doping. These results 





showed values of the minority carrier mobility that can 
only be explained with the incorporation of quantum 


mechanical scattering behavior at very high carrier 











concentrations. This appears to be the first direct 
measurement of diffusion lengths and minority carrier 
mobilities in this important material system [9]. Chapter 

develops the mathematical model underpinning the 














transport of minority carriers in the low dimensional 

















structures of interest. Chapter briefly describes the 





experimental apparatus and the technique used to extract 


the material properties, while Chapter IV explores the 

















theoretical limits of the model and demonstrates 
experimental evidence of those limits. In Chapter V the 
experimental evidence showing an increase in minority 








carrier mobility in heavily doped GaAs ~10°(cm*)is 





presented, and the results are discussed in the context of 


existing theoretical work. 





TH 


IS PAGE 








NTENT 





ONALLY LEFT BLANK 





II. TRANSPORT IMAGING IN THE TWO DIMENSIONAL LIMIT 


A. OVERVIEW 
SEM charge transport imaging combines two microscopes 
- a scanning electron microscope (SEM) to provide high 


resolution charge generation and an optical microscope to 





image the transport of charge. It can be performed in any 


material with a luminescent signature associated with 





charge recombination. In its basic operation, non- 
equilibrium minority carriers are injected into the 


luminescent semiconductor material by the SEM and the 











resulting radiative recombination is imaged through the 
optical microscope (OM). Analysis of the captured image 


allows quantitative, localized transport measurements. 





One application for this technique is as a means to 





perform contact-fr measurements of minority carrier 





diffusion lengths. This is a key parameter for many 


devices, including solar cells, photoconductors, and HBTs, 





as discussed in Chapter I. More conventional techniques 








for measuring diffusion lengths are generally limited by 
the need for contacts and the spatial averaging that occurs 


in macroscopic electrical measurements. Transport imaging 








can determine this important materials parameter directly 


from a single, Zero bias luminescent spot image, 





particularly for samples in the 2D (thin layer) limit. 


B. MODELING LIMITATIONS AND ASSUMPTIONS 





In specific application to the materials of interest 








from Chapter I we consider the case where a thin sample is 


doped and the charge generation rate is sufficiently low so 





that we are able to model the transport of minority 


carriers in an approximately constant distribution of 








majority carriers. For example, in doped samples of a lum 
active layer AlGaAs/GaAs heterostructure, with incident 
electrons of ~ 15 keV, this means an electron beam current 





of ~1x107° A or less for material doped at ~1x10'® om?. 
This approximation is made by assuming a highest value 
generation rate G of G ~ Eacc/Ei, where Eace is the incident 


electron energy and EB; is the energy required to produce an 








electron/hole pair. For energies in the ~ 5 - 40 kev 


range, one can approximate E; ~ 3E, for a bandgap of Ey 
[A215 The total carrier population An ~ Ap created then 


per electron is Gt. We approximate here a lifetime of Tl 





ns and a probe current of 1 nA, but the results can be 


scaled accordingly. In this example our ratio of resident 








majority carriers to minority carriers is on the order of 


100. Therefore, our low injection limit is valid. The 





generation volume radius for the electrons was approximated 


from the model of Kanaya-Okayama as ~ 1.5 um in GaAs at 30 


keV, with a hemispherical generation volume [13]. For more 





heavily doped materials, or shorter lifetime materials, one 
could use higher probe currents. Transport imaging can be 
performed outside these limits with more sophisticated 
modeling, but we will restrict ourselves to the low 


injection case for the analysis that follows. 


Cc. MATHEMATICAL MODEL DEVELOPMENT 





For cases where the diffusion length is comparable to 
or greater than our system resolution, diffusion of the 
minority carriers will broaden the luminescent spot. The 


extent of optical emission then becomes a function of 








minority carrier diffusion length and the diffusion length 


8 


can be directly extracted from the optical emission image. 





This approach cannot be easily applied to bulk/thick 
samples due to the generation volume created by incident 


electrons and the relatively weak dependence of the 





minority carrier distribution on diffusion length in 3D. 





However, since many new materials and devices utilize 





primarily thin films, (eg, heterostructures, quantum wells 
and specifically the base regions of HBTs) the range of 


applications for contact-free diffusion length measurements 





is large. 








In order to extract the diffusion length from the 
optical image, we model a steady state distribution of 
minority carriers created by a generation region of finite 
extent. The SEM beam, operated in a low injection 


configuration, is the source of the generation region, and 





our 2D assumption is based upon the relatively thin depth 





of the active region compared with its extent in the other 


two dimensions. 


25 keV e Beam 


e Penetration Radius 


GaAs Substrate 





Figure 3. E-Beam/Sample interaction schematic 








The sampl reaches steady state very quickly, and we 
will describe the distribution of the minority carriers in 
the optically active GaAs layer. The heterostructure has 
been modulation doped with Be (p-type), and the minority 


carriers are electrons. 


Beginning with the continuity equation for electrons 
in a p-type material: 


CLT A | en 


(1) aa a 
dt n n q n 


where G,is the generation rate UV ols 


An 
U, is the recombination rate = — for low injection. 
T 


n 


dn 
dt 





is the time rate of change of electrons per volume 


per second. 
J, is the current density vector 


An is the number of excess minority carriers available 


for recombination 


¢ 18. the Litetime tor electrons 





Defining the steady state current density for 


carriers: 





= = = Cc 
cm Ss 


HH, -ws- the mobility of electrons: in iGaks. and 


E is the externally applied electric field. 





10 


D, is the diffusion coefficient for electrons further 





related to the Diffusion Length by: 
(3) L=.Dt 
By combining the equations above we get: 


2 
(4) BG, 249. une En] 
dt T T 


n n 


By our assumption we are at steady state and therefore 


the time rate of change of the electron distribution is 





zero. Now, assume a constant E field in the x direction so 


that E=qEx, and Equation (4) becomes: 


2 
(4.5) 026 ap ey 
t dx T 


n n 


By multiplying through by A and making the 


substitution: S=yurE, where S is the drift length, we get 





Equation (5). 


= S 1 ord 
(5) Vintan, pea 





1. Generation Region 





Here we mak a digression to discuss the nature of 


Gt,, the steady state generation distribution created by 





the balance between the continuous SEM injection and 





recombination within the sample. After Donolato and 
Venturi we can define the distribution as a depth dosed 


Gaussian distribution [15]. 








(Z| = 
g(r, z;R)= 8o AAS 420°(R) 


6 
rn R 20° (z,R) 


1: 





Her th key featur of the distribution is the 








variance Oo being formed of two linearly independent 
factors: 
(7) o°(z,R) =o,’ +o0,"(z,R) 


where o, is the variance of the -béam and -o,. is: the spread 


of the primary and secondary electrons in the sample due to 





scattering. z and R are the depth coordinate and the 


primary electron rang respectively. R is a function of 





beam energy and the atomic number and density of the target 


material. 


Let. o, be the measure of the diameter of the circle 





within the beam that contains 50% of the total beam 
current, or d=beam diameter=1.67o0,. From empirical 


measurements we can assume that the lateral scattering 


3 
variance o,*0.15 [15]. R can be determined from ref [16] 





for a beam energy of 15keV in GaAs to be approximately 





xl5um, and =2.0um for a 25 keV beam energy. Though there 
is variation in the generation shape as shown in Figure 3 
as a function of the depth (z) we approximate the variance 
as constant for the generation region in our active layer. 
Assuming an average depth of R/2 the generation variance 


becomes: 


R 2 
(8) o (z,R)= 364°+0.1{ 2 





For a 15 keV beam energy and d=1.75um, o,,,=1.06um, and 


the radius within which 99% of the charge will be generated 
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define the limits of our 





is 2V20=3um. This 99% value will 


source region for numerical integration in later sections. 





We define our generation function for the source term 


using Equation (6): 





and normalize the output to 1 
(9) g(r’)=e 


2. Green’s Function Solution 
Equation (5) with the 





Returning to the differential 
Equation (9): 





inclusion of the source function, 





—n=— 


=) 8 
(10) Vint aan, — 7 7p 


By the use of an integration factor we can conduct a 





change of variables to eliminate the n, term thereby making 


Equation (10) into the Helmholtz equation. Substituting: 





(11) n(x y)=wixy,ye™ into (10) and combining terms we 


get: 
2 





ax Se ay ; 
(12) e (w.. tw, +2aw, ta’w) +e (w taw)——Fe Me ae 


Combining like terms and dividing through 


exponential: 





S , Sal) 15 4 
We PMG EW, POL |W Beda al ee 


in order to eliminate 


S 
Now we choose 2a+— =0 or a= 7 
L 2L 





the derivative term to obtain: 





S°+40 ee 
(13) W,, + Wy, ae Wes a 
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which is nothing more than the Helmholtz equation where the 


VS? +42 


Helmholtz operator is : V’+K’ and K=i oP 


Recognizing that the Green’s Function for the 





Helmholtz operator is the zeroth order Bessell Function of 


the second kind [16], 





(14) G(rsr')=5—K, (ker) 
V8? +42 


where k is the real part of K: k= an 


The general solution to a Green’s Function problem is 


Sx! =r? 


f 1) pee eee oh 
w(x, y) = [Gc re e*° dr 
0 


Returning to the solution for the electron distribution (n) 
= 
by the substitution: n(x, y)= w(x, ye“ or n(x, y)=w(x, ye 


Now, 


117, | vS?+42 ae ee 
(15) n(x, y) = ep. K, ———|r-r' 2L e2k eo dr’ 
2m Les OL 





By defining our limits of integration equal to the 


radius which contains 99% of the generation region’s 











minority carriers as described by the Gaussian distribution 


in Equation (8) our model for the minority carrier 





distribution including diffusion and drift in 2D becomes: 


220 2 2 S(x'-x) =r? 
{ K, VS°+4L0 ai 


2 2 
aF \r 2 e2a° dry’ 

















16 y= 
aes ae 


0 


The algorithm in Appendix A.5 is a Matlab coded version of 


this solution using a double quadrature numerical 
14 


integration scheme to calculate the distribution of 





carriers as they drift/diffuse from the finite area defined 


by the Gaussian generation function (9) and bounded by the 


circle of radius De. 
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IIIT. EXPERIMENTAL APPARATUS 


A. INTRODUCTION 


The minority carrier distribution in a semiconductor 





sample reveals information regarding the transport 








properties of the material itself. As shown in the 





previous chapter the diffusion length of the minority 








carriers determines the shape of that distribution. The 





novelty of the Transport Imaging technique is the 








extraction of the salient aspects of that distribution from 








an actual sample in a controlled and flexible manner. By 


combining the charge injection and high resolution electron 








imaging capabilities of a Scanning Electron Microscope 





(SEM) with the optical resolution of a Silicon Charge 








Coupled Device (CCD) Camera, accurate spatial 


representations of these distributions are captured and 














analyzed. A custom software solution allows for the 
analysis of these images and the fitting of the 
experimental data to the mathematical model’s predictions. 


B. APPARATUS DESCRIPTION 

The charge transport imaging instrument combines two 
independent microscopes - a JEOL SEM (See Table 1 for 
instrument specifications) for generating charge and an 


optical microscope (OM) for collecting and imaging the 





luminescenc mitted from the recombination process. Using 








a retractable arm, the OM is placed directly under the pol 
piece in the SEM, allowing the electron beam to pass 


through the center of the first optical collecting surface. 





The initial demonstration system modifies an OM attachment 





for the JEOL SEM that was originally designed to allow the 


LF 


fine adjustment of sample height required for wavelength 
dispersive X-ray spectroscopy (WDS). WDS provides higher 


energy resolution than traditional energy dispersive X-ray 








analysis (EDX) and requires sensitive control of sample 





height in order to maintain proper conditions for multiple 
diffraction angles. OMs designed for this purpose have 


short focal lengths and are normally used in conjunction 





with a lamp source and low sensitivity near-IR/visible 


imager. 


Trigger control 
= 










Optical 
Microscope 








Monochromator 


and PMT Image 


acquisition 
and 
display 











Retractable 
CL option Sample 





















































Continuous flow He 













































































Figure 4. Transport Imaging System Schematic 





A schematic of the system is shown in Figure 4. In 
addition to the charge transport imaging microscope, the 


instrument is equipped with standard CL capability using a 





Gatan (formerly Oxford Instruments) system with a parabolic 
mirror, ™“™ m monochromator and TE cooled GaAs PMT as the 
detector. Beam blanking capability exists for future time 


resolved work. Finally, the instrument has a liquid helium 





cooled stage, so that transport imaging and conventional CL 
can be performed over a temperature range from 300 to ~ 5 


K. The system uses continuous flow liquid He, with the 
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sample stage inside the SEM mounted in a cold finger 


configuration. 


Equipped with EDS capable of detecting Carbon and forming 
X-ray maps of composition; composition to within 0.1 wt% 


Integrated digital imaging system 
Noise reduction through frame averaging 


Image capture and export in electronic form (TIFF 
Low cost, medium quality thermal printouts 


High quality, medium cost Polaroid type 55 film containing 
both negative and positive 


Variable accelerating voltage; 200 to 40,000V 





Table 1. Table 1. JEOL 840A Specifications (From [11]) 





For transport imaging, the OM is used in a passive 


detection mode, detecting light emitted directly from the 





sample using a high sensitivity cooled Si CCD array camera. 
The current camera uses a 2184 x 1472 pixel array (15 mm x 
10 mm), with a pixel size of 6.8 x 6.8 uum? and can be used 


for transport imaging for wavelengths from ~ 350 to 1100 








nm. Initial image processing is performed using MicroCCD, a 





software program provided with the CCD camera. Although 


further image and data processing are often required for 








individual investigations, we benefit from excellent 





existing image acquisition and processing capabilities, 





often developed to support astronomical communities using 


Similar cameras for low light imaging. 
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The optical microscope insert is a basic two lens 





system (objective and eyepiece) modified to allow for 





passage of an incident electron beam. The considerations, 
as with any optical microscope, are resolution and 
magnification. Estimating the resolution for incoherent 


emission as 


0.614 
17 Dy os 
(17) a4 WA 


(where Ay is the spatial resolution, A is the wavelength and 








NA is the numerical apertur (set here at 0.95 max)), we 
find Ay = 0.56 um for A = 870 nm (e.g., room temperature 
emission from GaAs) and Ay = 0.22 um for A = 350 nm (e.g., 
emission from GaN). The current magnification of the 


optical system is ~ 20x, i.e., a5 x 5 um? area scanned by 





the e beam creates a 100 x 100 um? area on the CCD area. As 


mentioned, pixel dimensions are 6.8 um, so the resultant 





effective scale for the final image is ~0.4 um/pixel, 





comparable to the resolution limit for red/near IR light. 








In order to select photon emission from specific regions 
within the sample, appropriate combinations of bandpass 


filters are placed within the optical path for wavelength 











selection. The filters are used to eliminate, for example, 


substrate luminescence or to select the transport of 





interest in a multilayer sample. 


While the optical resolution limit is the fundamental 





mechanical limit of the luminescence collection, there 


exists a further analytical bound on the extraction of 





transport properties related to the data extraction method. 
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Cx DATA EXTRACTION 
As discussed previously, the transport property 


information is contained in the distribution of minority 





carriers at steady state. This distribution is directly 
linked to the resulting photon distribution as captured in 


an image by our device. Figure 5 shows one such image. 


‘Transport Imaging Luminescence Spot 


Pixel Number (Y-axis) 





80 100 120 
Pixel Number (X-axis) 


Figure 5. CCD image of Experimental Sample 


The data underlying this image is a 2x2 matrix whose 





indices correspond to pixel number. The values entered in 
each element (0-10, 000) of this matrix are the raw 
intensity of the photon emission collected by the CCD. In 


order to study the full extent of the distribution of 
minority carriers with greatest resolution, we extract line 
segments that cross the peak intensity point of the image. 
Though various methods may be employed to select these data 
sets, in this work that extraction was conducted via an 


algorithm written in MATLAB’ code (See appendix A.1- 








“imagedatamanipulator.m”). Once a line segment is 





extracted, it must be parameterized and fitted to the model 
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equation developed in Chapter ; This fitting can be 
accomplished via two methods with varying degrees of 


flexibility and resolution. 


1. Slope Analysis Estimation 
By assuming that the argument of the Bessel function 
is large compared with 1 we can assume the distribution 


approximates a negative exponential. 


=F 


7 a 
With r>>Lq : K( >e" 
‘d 


so that the slope of this line segment plotted on a semilog 


—1 
plot would be m= = c where 
L, MtkT 





m is the slope and r= x+y’, and all other terms are as 














defined in Chapter F Figure 6 shows a semilog plot of a 





line segment extracted from the image of Figure 5. 


LN(Iintensity) of Extracted Line Data 


LN(Normalized Intensity) 
& 


a AS = 05 0 0.5 1 15 2 
Radial distance from beam center (cm) x10" 


Figure 6. Semilog plot of extracted line of luminescence 
from Data Image of Figure 5. 
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Figure 7 shows the results of the slope analysis as 


calculated by the “Slope2.m” algorithm of Appendix A.3. 


Slope Analysis Estimate of Ld 
T 








© Data Points 
— L actual Ld=2.4484um 
— L La(+)=2.6237um 
— Ll d{-}=2.2731um 

© Data Points 
— Ractual Ld=2.4378um 
— RLd(+)=2.6139um 
— RLd(-)=2.2617um 











LN(Normalized Intensity) 


45+ 














ir ! i nt L | 
= 2 05 0 05 
Radial distance from beam center (cm) 3 


Figure 7. Slope Analysis for Experimental Sample 





The error bars here ar derived from the slope 
calculation of the distribution when the maximum possible 
mechanical error limits (104um) are assumed for the first 
and last points in the data sample. This error analysis is 


utilized instead of a standard deviation of data points 








from the linear regression due to its physical nature. It 
provides an intuitive link between the analytical 


assessment and the mechanical limits of our apparatus. 





The benefits of the slope analysis technique lie in 


its direct extraction of transport properties from an image 





with limited fitting or data manipulation. Depending on 





Signal to noise ratio and sample luminosity it can provide 


knowledge of a material’s diffusion length over a roughly 








3pm? area. This material property resolution is limited by 


the optical resolution of the system (pixel width) in the d 
23 





direction, equal to 04um, and the number of pixels required 





A 


to conduct the linear regression analysis in _ the r 
direction. Figure 8 shows a pictorial of resolution 
dependence on data sample size and _= selection. The 


resolution listed in the figure follows from Figures 4-6. 






Data Selection Region 
= Inormatized(!min)=- 1 0 


Data Line Segments 
g 5 Inormalized(!'max)=-01 





Luminescent Spot 


Figure 8. Slope Analysis resolution 





So in this case, the sample resolution of transport 





properties is averaged over an area of 0.4um x O6Oum or 
24um . 


As mentioned, the sample size and region selection is 
a function of the interplay between signal to noise ratio, 
error analysis and the large r limit. Lower error 


estimations require a larger number of data points, while 











for most samples, noise limitations drive our outer limit 
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below the ideal r/Lg>>1 limit. In work done by M. Talmadge 





of Fairfield University, Fairfield CT, it has been shown 





that when r>9Lqg the extracted -1l/m is within 95% of the 
actual Ig. Figure 9 shows the trend of predicted lg vs. 





actual Lg as a function of r/lLg[14]. 
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Figure 9. Slope Method Assumption Dependence on large 
Bessel Function argument 


However, when Lg is not known it is more difficult to 
determine this confidence factor. Work is currently being 
done to perfect a second derivative analysis to determine 
this confidence factor without a priori knowledge of the 
actual Lg. Additionally, when our samples are of low 


luminosity and have short diffusion lengths the collected 
25 








photon emission does not possess th required xtent to 


allow data selection within reasonable limits of r/L,g>>1. 





In this case a drift analysis is preferable however, it 








also possesses Similar limitations. 


2. Least Squares 2-Parameter Fit Analysis 
By applying an iterative least squares analysis (See 


algorithm Appendix A.4) one can fit the model prediction 








from Equation (16) to the full distribution of the 


extracted line of data and determine the diffusion length 








of the minority carriers and the radius of 99% charge 


generation. This analysis technique provides an additional 





understanding of the generation region dimension and is 





unencumbered by the limitations of the rfr/Lg>>1l limit. 








However, it increases the area of the sample used to 





extract a diffusion length ffectively reducing the 








resolution of the technigue from = 3 uum? to as large as 


~l6um> for the same sample data from Figure 5. An example 
of a completed fit is shown in Figure 10, where n is 


defined as the radius which encompasses 99% of the carrier 


generation or n=2N20 . 
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Normalized Experimental vs. Model Fit 





1 T T T T T T 


* Model Fit:(n=2.2um: Ld=2.2065um) Residue=0.19099 
— Data: S6 -PC=6e-10Amps 











Normalized Intensity 








Radial distance from beam center (cm) x10 





Figure 10. Least Squares Model Fit of Experimental Data 


While the area over which the parameters are 
determined is larger, this technique benefits from the 
possibility of greater accuracy in determination of those 
parameters. Moreover, when the material studied has a 
diffusion length greater than our optical resolution, with 
appropriately taken data and reasonable signal to noise 
ratios we expect to extract diffusion lengths and 
generations region radii accurate to within +0.lum. As 
mentioned before, this method is not limited by the 
necessity to take data far from the generation source, 


which is difficult for materials of low luminescence. 





However, ther are inherent limits to the materials and 





conditions which can be treated with this analysis. 
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IV. TRANSPORT IMAGING PREDICTIONS AND LIMITATIONS 


A. ERROR ESTIMATION FOR CURVE FITTING ALGORITHM 


Apart from noise and resolution limits, there are two 





primary sources of error in this technique; 1) the fit of 








the model data to th xperimental data, and 2) the 


assumptions that underlie our mathematical model. We can 





treat the quality of the model fit through the calculation 


of a root mean square error (RMSE). 


The residue listed in the legend in each figure is the 
sum of the least squares difference used to select the most 


appropriate parameter fit. From this residu on can 





calculate the RMSE of the fit, or 


d 2 M =1000 2 
ruse -( 4 where =residue= > (Lovet ~ Lana) and M is the 


i=l 


number of data points taken from the sample and used in the 








mode] calculation. This is a direct measure of the 








undetermined error of the fit to the distribution. Using 





this formulation we routinely achieve very favorable RMSEs 


(<10°). 


As was noted, our optical resolution is approximately 








0.4um, which results in only 100 data points taken over th 








40um interval shown in the figure. In order to smooth the 


distribution, we use a spine interpolation technique to 





increase our data set by a factor of 10. For relatively 


well behaved distributions, which these are, this technique 











has been shown to not alter the predictions of the model, 


yet allows a much higher confidence in the curve fit. 
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Becaus th model distribution is generated by a 








numerical integration calculation of an integral expression 





and the parameters are themselves arguments of non-linear 





functions, it is not readily apparent how this RMSE 





correlates to error bars in the parameters themselves. 
This analysis can be done and is explained in many non- 


linear least squares fitting texts, but it is more physical 





to vary diffusion length and generation radius by a small 




















amount and observe the resulting magnitude of the RMSE from 
our model fit. The next series of plots show this 
estimation. 

In order to establish a base line algorithm precision 
we first produce a model output for parameter values 








n=3.0um,L,=2.00um, then allow the least squares fitting 








routine to analyze that output and fit it with the 


appropriate parameters. When the integration step size is 





identical for the model produced output and the least 
squares fitting algorithm the residue is 0 as expected. 


This is shown in Figure 11. 
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Normalized Experimental vs. Model Fit 


Model Fit:(n=3um: Ld=2um) Residue=0 
Data: n=3L=2:PC=3e-3Amps 


Normalized Intensity 





2 -1.5 -1 -0.5 0 0.5 1 15 2 
Radial distance from beam center (cm) <1" 


Figure 11. Baseline calculation error in Model Fit 


algorithm (no variation of integration step size) 


A slightly more realistic assumption is that our 
numerical integration step size is on the order of 1000 


times smaller than that of the experimentally captured 





distribution. Figure 12 shows the effect of a difference 
in integration step size of a factor of 1000 between the 


model produced curve and the least squares fitted curve. 
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Normalized Experimental vs. Model Fit 


* Model Fit:(n=3um: Ld=2um) Residue=0.004692 
— Data: n=3L=2:PC=3e-8Amps 


Normalized Intensity 





2 -1.5 -1 -0.5 0 05 1 1.5 2 
Radial distance from beam center (cm) a6" 


Figure 12. Baseline calculation error in Model Fit 
algorithm (variation of integration step) 


The Residue of .004692 yields an RMSE of .0021 - an 
order of magnitude less than we typically see in our 
fitting to experimental data. We can, of course, reduce 
this step size but at the expense of integration time and 


given the very real effect of noise in our collected data, 








it is unrealistic to believe that we can achieve greater 


precision at room temperature through further reductions. 





In order to see the effect of parameter errors on the 





Residue we again employ the algorithm to produce an ideal 


model distribution and fit it with our least squares 





routine. The figures that follow show forced errors in the 


model fit of the same ideal data set and the impact on the 





size of the residue. We vary diffusion length (0.lum and 





0.2um) and generation region radius (0O.uwm and00lum) and 





calculate the RMSE in Table 2. It is noted that the RMSE 
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value is approximately the same for equal variations in 


either parameter. 


* Model Fit:(n=3um: Ld=2.1um) Residue=0.038476 
‘3e-BAmps 


Data: n=3L=2:PC= 





Figure 13. Error estimation for O.lum variation in Diffusion 


Length RMSE=6.2x10° 


* Model Fit:(n=3.1um: Ld=2um) Residue=0.034577 
Data: n=3L=2:PC=3e-8Amps 





Figure 14. Error estimation for O.lum variation in 


3 
Generation Region radius RMSE=9-8x10 
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Model Fit:(n=3.01um: Ld=2um) Residue=0.0064142 
Data: n=3L=2:PC=3e-8Amps 





Figure 15. Error estimation for 00lwm variation in 


Generation Region radius RMSE=2.5x10~ 





Another reference point for observed errors is shown 


Figure 16. 


= Model Fit:(n=3um: Ld=2.2um) Residue=0.11098 
Data: n=3L=2:PC=3e-8Amps 





Figure 16. Error estimation for 0.2um variation in 


Generation Region radius RMSE=10.5x10° 
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in 
































Residue RMSE L or n Variation 
.0064 2.5x10° +0.01um 
.0384 6.2x10~° +0. 14m 
.1110 10.5x10° +0.2 um 
.4763 21.8x10° +0.5 um 
1.6853 41.1x10~ +1.04m 
Table 2. Tabulated Error Estimates for Curve Fits 
B. LIMITS OF MODEL ASSUMPTIONS 


A more fundamental error resides in the boundaries of 
where our model assumptions break down, or where other 
aspects of transport begin to play a more dominant role. 


There is much interesting science in this aspect of the 





analysis, and in fact, Chapter V will focus on the 
interplay of one such phenomenon, photon recycling, which 


at high doping levels begins to affect the luminescence 





distribution on a scale that demands special treatment. 


1. Low Injection Assumption 


An important limitation in our modeling is our 

















assumption of low injection. As described in Chapter Bis 





for these samples we are restricted to probe currents equal 





to or below I1xl0°A. Above this level of excitation we 


Significantly alter the distribution of majority carriers 








in the vicinity of the generation region and the 


recombination is no longer appropriately described as 





proportional to the density of minority carriers alone. In 


order to probe this limit and to compare the slope analysis 





predictions with the model fit technique, we will observe a 
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series of data taken from the same spatial location on an 





AlGaAs/GaAs heterostructure kept at a temperature of 4.7 K, 


and constant beam energy of 25 keV while the probe current 








was varied from 6x10°’-6x10°A. Previous work within our lab 
has reported th ffect of increasing SEM probe current on 
the size of the luminescent spot [17]. Here a series of 


images is presented corroborating this work and quantifying 
the increase in the standard deviation of the generation 
distribution. Figure Le f shows a schematic of the 
heterostructure as designed by Tom Boone at Hitachi Labs. 


Gao 6Alo.4As: 0.2um; p-5x10'® cm? 
electron confinement 


Grading: 500 A interface recombination 


Gao ¢Alo.4As: 0.1um PL active region ~870nm 
Na-5x10'® cm @ 


Grading: 500 A interface recombination 


Gao.¢Alo.4As: 0.2um; p-5x10'® em? 
electron confinement 





Figure 17. AlGaAs/GaAs Heterostructure design from Tom 
Boone Doctoral Dissertation [22] 














The study of these materials at low temperatures allows a 


greater signal to noise ratio, which enables greater 





accuracy of the model fit. Previous work in our lab has 





shown that the minority carrier diffusion length in these 


materials is independent of sample temperature, and 








therefor W can us thes measurements to establish a 
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baseline of accuracy between the two techniques that should 






































translate to higher temperatures [23]. Table 3 compiles the 
salient results from the analysis of this sample at 4.7 K. 
Probe La by BTePe ha by Model 2/20 TT bore RMSE of Model 
Current | Analysis (pm) | Fit (um) ; Fit 
Region 
6x10°7A | 1O+8 3.9+0.5 3.2+0.5 25x10 
6x10 'A | 2.9+0.2 3.7202 3.8+0.2 12x10~ 
6x10 °A | 2.9+0.2 4.1+0.1 4.2+0.1 9x10" 
6x10°A | 3.4+0.2 4.4+0.3 4.5+0.3 15x10” 
6x10°A | 3.440.2 4.7+0.3 4.9+0.3 15x10° 
Table 3. Measurement Results for O.lum active layer, Boone 





Heterostructure #9 


2. Slope Analysis Limitations and the Low Injection 
Limit 


As can be seen from the tabulated values, increasing 





probe current tends to increase th ffective radius of the 


generation region, as expected and reported previously in 














ref [17]. The diffusion length as measured by both 


techniques is relatively constant as a function of probe 





current below 6x10 °A in accordance with our model 





assumptions. Also evident is the disparity between the 
slope analysis method and the model fit. This difference 


is expected and is related to the degree to which the slope 








analysis limiting assumption of large Bessel function 


argument is valid. That is, the slope analysis predictions 





assume that we are in a regime where the Bessel function 
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r . 
argument —>>l. In this case, measurements were taken over 
‘d 


a distance of 7-14um from the center point. Assuming an La 





actual of 4um we can estimate the degr of disparity that 





should result by consulting Figure 9. Entering the X-axis 
with a value of 7/4or 2.5 we extract a Talmadge factor of 


0.75, or we would expect that the slope analysis method 





would predict a value within 75% of the actual. This 
corresponds well to the ratio of the model fit prediction 





to slope analysis prediction 3um/4um —- or 0.75. The lowest 
probe current shows the limits due to noise in this 


analysis. As expected, the slope analysis method will be 








impacted more greatly by poor signal to noise ratios 





because of its higher spatial resolution. 


The figures that follow show the comparison between 





slope analysis plots and model fit plots. Tet is 





instructive to observe which portions of the model fits 





begin to deviate from th xperimental data for higher 
probe currents. The trend away from low injection can be 
tracked by observing the deviation in the “shoulder” 


regions of the distributions as the probe current 


increases, (Figures 20,21, and 22.) 
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Figure 18. Slope and Model fit analysis plots for 6x10’A 
probe current (pertinent data tabulated in Table 3) 





Figure 19. Slope and Model fit analysis plots for 6xl0"'A 
probe current (pertinent data tabulated in Table 3) 


39 


Mote Feed Dae Lind Rt Fete THE 
— tae SF Poke these 





Figure 20. Slope and Model fit analysis plots for 6xl0°A 
probe current (pertinent data tabulated in Table 3) 
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Figure 21. Slope and Model fit analysis plots for 6xl0°A 
probe current (pertinent data tabulated in Table 3) 
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Figure 22. Slope and Model fit analysis plots for 6xl0°A 
probe current (pertinent data tabulated in Table 3) 


From Table 3 we can see that the balance between 


signal to noise, and model limitations place the best probe 


current for accurate measurement at 6xl0"°A for this sample. 





3. Small Diffusion Length Limitations and the Role of 
the Generation Distribution 


A more substantive measurement limitation for the 











heavily doped materials discussed in Chapter I is the 
relatively small diffusion lengths that accompany such 
large concentrations of acceptor dopants. While the 





literature predicts an increase in the minority electron 


mobility in GaAs doped with Be starting at 5xl0“cm”, [23] 


the lifetime continues to trend downward at a rate which 


overpowers the increase in mobility and causes diffusion 


lengths to continue to decrease. Just beyond this 


concentration the diffusion length drops below lum, and we 





approach another limit of our technique. This limit is 
In order to see 





directly related to the generation region. 
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how this limit arises, we again study our solution to the 
transport equation for a Gaussian generation distribution. 
Ly) =r? 
2/20 [524 47? S(x'-x) 


1 ’ 2 o ' 
(16) MOD) = ee [ K, P Ir-r'| IP” pio” dy 


0 











Here the two terms which contribute to a zero E field 


(S=0) diffusion ar th Bessel function and the Gaussian 





source function. As previously discussed, our use of the 


Bessel Equation comes from well-established differential 








equation theory for solving diffusion equations in 2D 





carrier transport and other disciplines governed by the 


Helmholtz Equation. The assumption of the Gaussian 





distribution to represent the carrier generation region 


within the sample is based upon the statistical 








interpretation of the electron-electron scattering and in 
limited cases is backed by empirical evidence 


[18], [20], [21]. 





More recently, work in our lab has shown that these 








distributions may not all fit the same mathematical 


dependence [19]. By assuming a standard Gaussian 














distribution, we may be neglecting effects of small 
deviations due to sample geometry, beam inhomogeneities, 


and possibly other effects that govern the sub-micrometer 





scale granularity of the minority carrier distribution. 


These inaccuracies in our model will become more prevalent 








when materials of small diffusion length are studied. 





Moreover, it is anticipated that when materials which have 


diffusion lengths on the order of our optical resolution 











are studied, it will become difficult to observe th £LECE 





of diffusion on the distribution. At this limit we may say 
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that we are indeed observing the generation region itself. 











As we approach this limit, the accuracy of our 
representation of the interaction region will become 
increasingly important. Any deviation of the actual 
distribution from our Gaussian model will be reflected in 


some mixture of parameter adjustments, which will 





unrealistically be portrayed as diffusion length or sigma 


variation by the fitting algorithm. 





In order to make a quantitative assessment of this 





limit we will demonstrate a limiting case. By producing a 
distribution data set defined as a pure Gaussian (n=3.0um) 


and allowing the fitting algorithm to fit Equation (16) to 








it, we may see what a material with no diffusion and a 


perfect Gaussian generation distribution might look like. 





Figure 23 shows the slope analysis and model fit for this 
case and demonstrates that both methods inaccurately 


predict a diffusion length of 0O.3um. As predicted, the 





model fit compromises the generation region radius of the 


distribution from its known value of 3.0um in order to fit 





the data with a diffusion length that allows for the 


smallest residue permissible. 





Slope Analysis Estimate of Ld 
T 








Normalized Experimental vs. Model Fit 








© Data Points 
— Lactual Ld=0.41058um 
— LLdj+)=0.5891um 
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© Data Points 
— Ractual Ld=0.41058um 
— RLd(+)=0.5891um 
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Figure 23. Pure Gaussian Distribution Model Fit 
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In other limiting cases where the modeled data set was 





produced with very small but non-zero diffusion lengths, 





the algorithm’s accuracy was directly proportional to the 


step size, (analogous to pixel size) used in the creation 





of the data. This is consistent with our prediction that 





our CCD pixel size results in an effective lower limit of 











discernability 0.4um for either parameter value. However, 


because the generation distribution radius does not 








approach this lower value, it is effectively only a limit 


for our determination of diffusion lengths. Compounding 





this lower limit is any inaccuracy in our assumption of 





generation region form, which will tend to distort both o 











and lg. As we approach the regime where th form of the 
generation distribution contributes more than does 
diffusion to the shape of the overall distribution we 
expect this limitation to have a larger and larger effect. 





Moreover, the 0.4um error suggested above is only accurate 
when you assume a perfect Gaussian generation region. Any 


deviation of the generation region from the ideal will tend 





to increase our baseline error. 





In reference [19] Luber discusses the use of this same 








Transport Imaging technique aS a means to more accurately 








determine the interaction region distribution for materials 








of interest. While a full quantitative method has not been 


developed, it is seen as a key step toward the study of 





very small diffusion length materials (L,<1.0um). 
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V. STUDY OF HEAVILY DOPED HETEROSTRUCTURES 


A. MOTIVATION 





As discussed in Chapter I, faster switching 
transistors are OF prime importance to military 


applications. The devices that are currently under 





development to handle this task, in commercial applications 





as well as military, are HBTs. The key device parameter 





for increasing speed and efficiency is the base layer 





transit time. Many engineering design approaches are used 





to decrease the time it takes electrons to flow across the 





base layer, but the efficient use of these techniques is 
dependent upon the accurate knowledge of the transport 


properties - minority electron diffusion length, lifetime, 








and mobility. 


A mathematical description of the base transit time 
reveals the importance of low dimension construction as 


well. 


(18) Tt, == 


where ft, is the base transit time, Wis the base width, and, 


(19) p, =HlT/ 


is the minority carrier diffusion coefficient as defined by 





the #instein equation where jf, 1s the minority carrier 


mobility, e the charge on an electron, k Boltzman’s 


constant, and T is the temperature in Kelvins [1]. 





The appearance of the base width as a squared term 


dominates the trend of transit time, but with decreased 
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base width comes the problem of emitter-collector current 
leakage. Here, increased dopant concentrations aid the 


reduction of the base width by providing an impediment to 














this current leakage, however, classical analysis predicts 





that increased doping has a negative effect on 7, through 





its reduction of mw, and therefore D,[1]. Classically, one 





would predict a practical limit to the concentration of 


dopants that can be used as the competing effects of 


reduced D, and decreased W interact. However, a more 


detailed analysis reveals a more complex picture. 


B. QUANTUM MECHANICAL PREDICTIONS 
An observation of nonconventional electron current 


density in an GaAs/AlGaAs N-p-n HBT with Be base doping of 


6x10°cm? led Lyon and Casey to believe that some other 








transport mechanism was at play. They observed a 
collector-emitter current density that exceeded 
conventional predictions by four times [29]. More 





recently, material growth techniques have improved and base 


layers are being produced with graded doping schemes in the 








low 10%°cm? [5]. At these levels the assumptions of the 
Boltzman distribution and classical carrier scattering 


descriptions may not be sufficient to explain carrier 





transport. In work done by Bennett and Lowney employing a 





first principles quantum mechanical analysis of scattering 





mechanisms in heavily doped GaAs, it is predicted that a 


local minima exists for electron mobility in the 5xl0°cm> 


regime. This analysis includes all the important 





scattering mechanisms for the low-field mobilities: 





acoustic phonon, polar optic phonon, piezoelectric, ionized 
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impurity, carrier-carrier, alloy, and plasmon scattering. 





The upturn in the mobility results from the dependence of 





these scattering mechanisms on the dopant and carrier 


density. As the dopant density increases the average 





distance between holes decreases. This screening radius 


then determines the upper frequency that may be supported 








for the vibrational modes set up in the plasma of majority 
carrier holes, plasmon cutoff frequency (PCF). As the PCF 


increases, the scattering interaction probability between 











minority carrier electrons and plasmons decreases. 
Additionally, as the fr hole concentration increases the 
lower energy bands fill, and Pauli Exclusion Principle 
screening becomes important. The number of majority 





carrier/minority carrier scattering events is reduced 





because the holes are precluded from changing their energy 


level and therefore can not interact [23]. 


These results have been difficult to reinforce 





experimentally and limited direct videnc xists to 





support them [25-28]. A method to observe this effect in a 


non-contact manner which requires little sample preparation 





and is non-invasive would add to existing device 








diagnostics techniques. Transport Imaging provides such a 


solution. 


Cc. EXPERIMENTAL RESULTS 





In order to observe this increasing mobility trend a 





series of Be doped heterostructures was studied. Figure 24 





shows the design of the studied structures. 
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Gao. ¢Alo.4As: 0.2um; p-5x10"8 em 
electron confinement 


Grading: 500 A interface recombination 


Gao ¢Alo.4As: 1.0um; PL active region 
Na-3x10'® -1x10° cm?  =peakemission ~870nm 


Grading: 500 A interface recombination 


Gao.¢Alo.4As: 0.2um; p-5x10'® em? 
electron confinement 





Figure 24. Heterostructure design as grown by Tom Boone, 
(From Ref. [22]) 





The Transport Imaging technique was applied to seven 
samples of differing active layer doping and the diffusion 


length of the samples was extracted using both the slope 





analysis method and the two parameter fit. Data were taken 


at multiple locations on each sample, on different days and 





under varying beam energy and probe currents. Though some 


variation of parameter value with location was noted, 





overall the samples can be considered to be very 


homogeneous, and these results to be representative of the 





average properties. In order to test the predictions of 








Dr. Bennett, we required lifetime (t) values with which we 


could extract mobility (uw) values from our diffusion length 





measurements through the relationship of Equations (3) and 





(18). We coupled independent measurements of the sample 








lifetimes with values provided by the grower, Dr. Boone. In 





both cases th measurements were made by time resolved 





photo-luminescenc techniques. Table 4 tabulates the 
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initial measured values 


marked with an asterisk wer 


for these 


measured 





samples. 





Th 


lifetimes 


in 2003-2004 by Yale 





University and were not corroborated by our independent and 


more accurate TRPL confirmation. 


Yale repor 





ted lifetimes and our TRPL measurements for 








Comparisons between th 





the 


















































other samples showed a 10% overestimation in the samples of 
with Na>x10?? cm®. 
Sample Doping[cm™°] t [ps] La [um] fon" 
OF a eS 
Ves 
A2 2.75x10'° 2050 * 3.64.1 2500+140 
B7 3.75x10'8 900 2.34.1 2350+ 200 
Gg 5x10'® +3x10!8 4800 4.1+.1 1333+ 68 
D6 35x10" 95 1.64.2 10800 + 2800 
E3 4.0x10'° 140 * 1.84.15 9200 +1600 
F4 6x10” 116* |e aes | 12400 +1300 
G8 1.0x107° +.1x10”° 11 1.74.2 101000 + 25300 
Table 4. Initial results of AlGaAs/GaAs heterostructure 
study 
1. Initial Observations 


The result for sample C9 
completely consistent with 


is also the sample tested in Chapter 


active 


samples. 


trend of decreasing mobility 


a concentration of 5.0x10%cm* 


2 
1000. 
V 


‘S 


layer 


dimension 


10 times 


Drs 





Bennet 


is the only sample that 


























Other reported data for mobilities 


thinner 


than 


in 


is 





t’s predictions. It 


, and possesses an 


the other 





Samples B2 and C7 are consistent with Bennett’s 


toward the inflection point at 


though offset by approximately 
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this range 





of doping concentrations also show elevated values compared 


with Bennett’s predictions [24]. However, the mobilities 





for samples D-G are unrealistic, even if the trend of 





increasing mobility is evident. 


The consistent positive offset of these values 
suggests a systemic error or effect that is operational in 
the: S10" "em - samples. We propose two reasons for these 


offsets and apply appropriate offsets to account for them. 


2. Generation Region Discrepancies 


As mentioned in paragraph B, we expect that as we 





approach the regime where the generation region contributes 
more and more to the shape of the extracted curve, we will 
be subject to limitations due to inaccurate assumptions 
about the generation region. By studying the curve fits in 
the higher doped samples we can gain some intuition about 
where this limit may be. The figures that follow (Fig. 25- 


28) are the best residue curve fit achieved for each of the 





samples. 











Normalized Experimental vs. Model Fit Normalized Experimental vs. Model Fit 


* Model Fit:(n=3.5um: Ld=3.6um) * Model Fit:(n=3.7um: Ld=2.28um) 
09 Residue=0.035 oot Residue=0.055 
i — Data: $20 -PC=6e-10Amps 3 — Data: $72 -PC=6e-10Amps 
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Normalized Intensity 
° 
on 
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Figure 25. Samples A2 and B7 best fit 2-parameter fit 
extractions 


50 


[+ Mosel Fit (m=2 Sum Lael oa ide, 
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Figure 26. Samples D6 and E3 best fit 2-parameter fit 
extractions 


* Model Fit:(n=2um: Ld=1.9um) 
Residue=0.057 
— Data: S40 -PC=6e-10Amps 





Figure 27. Sample F4 best fit 2-parameter fit extraction 


* Model Fit:(n=0.8um: Ld=1.89um) F * Model Fit:(n=2.4um: Ld=1um) 
Residue=0.12 i Residue=2.5 


— Data: $8 _PC=e-10Amps — Data: $8 -PC=6e-10Amps 





Figure 28. Sample G8 best fit (a) algorithm run, (b) 
Assumed reasonable generation region with algorithm fitted 
diffusion length 


eal 


The generation region radius has been shown to 
increase with probe current and with beam energy, but here 


we s a decrease in the generation region as a function of 





dopant concentration. Also noted is the tendency of the 








model fit to depart from the data set in the region of the 





shoulders of the curves near the base of the distribution. 
This effectively causes an overestimation of the diffusion 


lengths as predicted in Section B. Finally, we can surmise 











that in samples D-G we are in the realm where the diffusion 


lengths are on the order of the generation region radius, 





or less, and approaching a fundamental limit of our 





assumptions. If we assume that the produced error is on 


the order of that demonstrated with the pure Gaussian from 





Chapter IV, Section C. we would expect an overestimation of 


the diffusion length by 0.4um. 


3. Photon Recycling (PR) 





Another important and well-documented effect that must 


be considered is that of Photon Recycling. The literature 





is replete with documentation of this phenomenon that 


affects diffusion coefficients and observed lifetimes in 





bulk GaAs that begins to act in this doping regime [30-32]. 
The effect is treated in different manners, but 
consistently results in correction terms being used to 


adjust the observed diffusion coefficient and total 














lifetime. Renaud treats the effect as an addition to the 
generation function in the continuity equation (our 
Equation (10)). He defines the photon recycling generation 
function: [30] 
a; i U Ud Uy 
(20) Gyp => | KO x) An’, dx 
2% 


52 


This represents the excitation in the sample with thickness 
w and average absorption coefficient a, taking into account 


the spectral density of the light. Tr is the radiative 


lifetime related to the lifetime we measure with TRPL by 














1 1 1 
(21) —=—+4+— 
a a 
where Tae is the non-radiative lifetime. The real 


perturbation comes from the calculation of K(x,x’), which 


is related as a series of exponential integral functions 








[30]. The minority carrier distribution is then expanded 
in a series expansion over the photon recycling source 
region and the continuity equation is now adjusted with 


each term possessing a PR pertubation factor Ty, where: 


1 ae 
(22) jes } K(x,x+u)u"du 
nit 


ro-x 


Because K is principally a function of exponential integral 





functions and converges quickly to zero with increasing n, 





they can be represented by the spatial average value Ey 
Appendix B lists the first two non-zero terms of this 
Series: (T,)and (T)) in their full mathematical form, as well 


as the exponential integral function. The resulting 





continuity equation is a modification of our Equation 


(4.5): 





2 
(23) 0=6,-[+-(n)} nett +n) 
T 


Renaud demonstrates good agreement between his 


corrective terms and behavior of GaAs LEDs and photovoltaic 
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cells, and calls for more study on small thickness samples. 


Badescu states that while photon recycling is most apparent 





in bulk samples, there is a more pronounced effect in 


1 
samples where the absorption length L,=— exceeds the 
a 


diffusion length, even for thinner samples (<lum) [31]. 





If we add a column to Table 4 and populate it with the 
absorption coefficient a (taken from [33] ) and the 
absorption length for each sample we see a correlation 


between the departure from predicted values of mobility and 





the breakpoint where absorption length exceeds diffusion 











length. 
Sample | Doping[cm ’] t [ps] Tig oo t ene 
[wm] Veg 
bum [pm] [pum ] 
A2 2.75x10'° 2050 * | 3.6+.1 5000 2.0 2500 +140 
B7 3.75x10" 900 4500 2.2 2350 +200 





C9 5x10 +3x10'" | 4800 4000 | 2.5 1333468 





3500 | 2.85 | 10800+2800 





D6 3.5x10" 95 











F4 6x10" 116* 3000 | 3.3 | 12400+1300 


2.34.1 
4.1+.1 
1.64.2 

E3 4.0x10" 140* | 1.84.15] 3400 | 294 | 9200+1600 
194.1 
ee 


G8 1.0x10” +.1x10” 11 2300 | 4.3 | 101000425300 























Table 5. HS data table with absorption length comparison 





Cc. CORRECTIONS FOR OPERATIVE EFFECTS AND DISCUSSED 
LIMITATIONS 


The photon recycling effect is dependent upon the 





number of photons generated, the rate at which they 
reabsorb, but also on the rate at which they can escape the 


active layer befor creating additional electron-hole 














pairs. The first dependencies we have previously 
described, but now we must look at the index of refraction 


of our samples and the corresponding critical angle of 





total internal reflection. 








In 1976 Asbeck reported the critical angle for 


GaAs/AlGaAs interfaces as a function of various Al 





concentrations [35]. Interpolating from his graphs and 





confirming with Snell’s Law, we arrive at a critical angle 





for the active GaAs layer of 69.6° The index of refraction 


varies with dopant density as_ well, but becaus th 





differences are small between AlGaAs and GaAs we can use 


the value for GaAs as 3.59 and for 40% Al concentration in 














AlGaAs n=3.36. Using these values to calculate the 


Reflectance; [36] 


l-n) 
= a (i=4) 
+n 
NaiGaAs = 3.36 
399 





For n= =.935 and therefore R=0.1%, or when light 





"GaAs 





strikes the interface at an angle less than the critical 











angle, 99.9% will transmit through to the AlGaAs layer. 
Also required is the radiative lifetime. From [32] we can 
define: 
1 
(24) a Ls] 
210s N, 


Now to calculate the correction factors and apply them to 





















































our experimental results for the samples in Tables 3 and 4 
we employ the MATHCAD routine of Appendix B.1. The 
detailed calculation sheets are in Appendix B.x and the 
overall results are tabulated below in Table 6. 
Sample Doping Tt [ps] La** cm Tpr Lapr | com? 
pb | — Her | —— 
[om73] [um] V-s os [pm] [V-s 
A2 2.75x10'® 2050 * | 3.6+.1 2500 +140 NC NC NC 
B7 3.75x10" 900 its am 2350+ 200 NC NC NC 
C9 5x10'8 +3x10'° 4800 4.14.1 1333+ 68 NC NC NC 
D6 35x10" 95 12 6000 + 2000 69 10a. 6000 + 2000 
E3 4.0x10'° 140* | 1.44.15 5600 + 1300 85* | 1.14.15 5500 + 1600 
F4 6x10!” 116* | Ds see | 7800 + 1000 66* | 1.0+.1 7600 + 1300 
G8 1.0x107 +.1x10” 11 1.0+.2 | 59000+ 20000 9 0.9+.2 | 36000 +15300 





























Table 6. Tabulated Parameters corrected for Generation 
Region error(**) and Photon Recycling overestimation 
D. CONCLUSIONS 


The effect of the Generation Region (GR) error plays a 


much stronger role on the calculated mobility values than 











does Photon Recycling (PR) because of the simultaneous 
effect PR has on diffusion length and lifetime. The PR 
effect is seen to grow as a function of doping. 


The assumption of a 0.4um error for GR is an estimate 








that needs more refinement, through the development of an 


analytical assessment of generation region definition and 





its inclusion in the numerical integration algorithm of 
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ity (j2./ p20) 





normalized electron mobi 





Appendix A. Transport Imaging provides an appropriate 


mechanism for this analysis and should be pursued. 


The mobility values from the final corrected column 








are plotted against Dr. Bennett’s predictled results with 


appropriate error bars in Figure 29. 


Minority Electron Mobility ;€Table 2.1) 











6 
~ 
4 
3 
2 
1 
QO GaAs exp. data, Harman, et al, 
0 Appl. Phys. Lett. 83, 537, 1993. 
16 17 18 19 20 21 
° 3 
log (acceptor density - cm’) 
rsd = 0.018 polynomial/transformation 
Figure 1 
Figure 29. Final corrected Transport Imaging mobility 
values reported (After Bennett [9]) 
The local minimum is clearly demonstrated though 
absolute magnitude agreement is not. A new round of 


experiments is planned to test the magnitude relationships 
57 
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ce 





of the wt product through th measurement of Larize by 





studying the distributions as a result of an applied DC 


bias. In this manner a full distribution fit should escape 








the limitations resulting from generation region error, 





though the optical resolution limitation (0.4 mum) may still 





be operative in the samples of heaviest concentrations. 


From the results demonstrated it can be assumed that 








of the limitations and constraints inherent in Transport 








Imaging the assumption of a generation region distribution 





has the largest impact for measurement of low diffusion 





length materials. It appears that experimental results can 


be assumed valid so long as the diffusion length measured 





is on the order of the generation region radius (as in 








Samples A9-C2), that the signal to noise ratio is favorable 
(as in all data samples shown herein), and that the 


diffusion lengths measured are greater than the optical 





resolution of the system (0.4um in these data samples). 





Several methods may be useful to overcom thes 
constraints and are being studied in our laboratory. They 


include time resolved techniques reminiscent of the Haynes 








Shockley xperiment, but maintaining the spatial 





information of the light emission to great resolution, AC 





drift techniques attempting to generate resonance responses 





between transport properties and the applied electric 
force, and observation of the effects of magnetic fields on 
the flow of the charge carriers at the sub micrometer 


scale. 
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function varargout = newfirsttry(varargin) 


$29 August - changed mkdir location for gui directory creation to create 

within diffusive data analysis\analyzed data. 

Will change parameter matching to Ld versus Mobility 

$19 September - updated version on laptop with amended version from desktop 

%12 October - Cleaned up comments and reorganized order of functions 

$16 October - Reset all flag values in each main analysis functions (aflag,gflag, 
flag) all tested for empty variable status and set to '0' if so. 


3p 


NEWFIRSTTRY Application M-file for newfirsttry.fig 
NEWFIRSTTRY, by itself, creates a new NEWFIRSTTRY or raises the existin 
singleton’. 


% H = NEWFIRSTTRY returns the handle to a new NEWFIRSTTRY or the handle to 


the existing singleton’. 


NEWFIRSTTRY ('CALLBACK', hObject, eventData,handles,...) calls the local 


function named CALLBACK in NEWFIRSTTRY.M with the given input arguments. 
NEWFIRSTTRY('Property', 'Value',...) creates a new NEWFIRSTTRY or raises the 
existing singleton*. Starting from the left, property value pairs are 
applied to the GUI before newfirsttry OpeningFunction gets called. An 
% unrecognized property name or invalid value makes property application 

stop. All inputs are passed to newfirsttry OpeningFcn via varargin. 
*See GUI Options - GUI allows only one instance to run (singleton). 

See also: GUIDE, GUIDATA, GUIHANDLES 

Edit the above text to modify the response to help newfirsttry 

Last Modified by GUIDE v2.5 30-Mar-2005 12:09:12 


% Begin initialization code - DO NOT EDIT 


gui_Singleton = 1; 


gui_State = struct('gui_Name', mfilename, 
‘gui Singleton', gui Singleton, 
‘gui OpeningFcn', @newfirsttry OpeningFcn, 
‘gui_OutputFcn', @newfirsttry OutputFcn, 
'gui_LayoutFen', (l, 
"gui_Callback', El}? 


if nargin & isstr(varargin{1}) 
gui_State.gui_Callback = str2func(varargin{1}); 
end 


if nargout 
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fit? 
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49 varargout{l:nargout} = gui_mainfcn(gui_State, varargin{:}); 

50 else 

51 gui_mainfcn(gui_State, varargin{:}); 

52 end 


53 % End initialization code - DO NOT EDIT 

54 %Initialize global variables 

55 global aflag gflag foldername firsttime fitflag; 
56 %aflag=0; 

57 %gflag=0; 

58 %fitflag=0; 

59 %firsttime=0; 

60 filetoload=[]; 








61 
62 % 4 
63 % --- Executes just before newfirsttry is made visible. 


64 function newfirsttry_OpeningFcn(hObject, eventdata, handles, varargin) 
65 % This function has no output args, see OutputFen. 
66 % hObject handle to figure 


67 % eventdata reserved - to be defined in a future version of MATLAB 
68 % handles structure with handles and user data (see GUIDATA) 

69 % varargin command line arguments to newfirsttry (see VARARGIN) 
70 


71 = % Choose default command line output for newfirsttry 
72 handles.output = hObject; 

a3 

74 % Update handles structure 

75 guidata (hObject, handles); 


76 

77 if nargin == 3, 

78 initial_dir = pwd; 

79 elseif nargin > 4 

80 if strempi(varargin{1},'dir') 

81 if exist (varargin{2},'dir') 

82 initial_dir = varargin{2}; 

83 else 

84 errordlg('Input argument must be a valid directory','Input Argument Error” 
ae 

85 return 

86 end 

87 else 

88 errordlg('Unrecognized input argument', 'Input Argument Error!'); 

89 return; 

90 end 

91 end 


92 % Populate the listbox 
93 load_listbox(initial_dir,handles) 
94 % Return figure handle as first output argument 


95 
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96 % UIWAIT makes newfirsttry wait for user response (see UIRESUME) 
97 % uiwait (handles.newfirsttry); 


98 
I a Pg ee ce oe | 
100 &% --- Outputs from this function are returned to the command line. 


101 function varargout = newfirsttry_OutputFen (hObject, eventdata, handles) 


102 % varargout cell array for returning output args (see VARARGOUT); 
103 % hObject handle to figure 

104 % eventdata reserved - to be defined in a future version of MATLAB 
105 % handles structure with handles and user data (see GUIDATA) 
106 

107 % Get default command line output from handles structure 

108 varargout{1} = handles.output; 

109 

110 ¢$ 





111 % Callback for list box - open .fig with guide, otherwise use open 
112 % ------------------------------------------------------------ 
113 function varargout = listboxl_Callback(h, eventdata, handles) 


114 % hObject handle to listboxl (see GCBO) 

115 % eventdata reserved - to be defined in a future version of MATLAB 

116 % handles structure with handles and user data (see GUIDATA) 

117 

118 % Hints: contents = get (hObject, 'String') returns listboxl contents as cell array 
119 §& contents({get (hObject, 'Value')} returns selected item from listboxl 

120 


121 get(handles.newfirsttry, 'SelectionType'); 
122 global filetoload tiftoload ; 


123 

124 if strcemp(get (handles.newfirsttry, 'SelectionType'), 'open') 
125 index_selected = get (handles.listbox1, 'Value'); 

126 file list = get(handles.listbox1, 'String'); 

P27 filename = file list{index_selected); 

128 if handles.is dir (handles.sorted_index(index_selected) ) 
129 cd (filename) 

130 load_listbox (pwd, handles) 

131 else 

132 (path, name,ext,ver] = fileparts (filename); 

133 switch ext 

134 case '.fig' 

135 guide (filename) 

136 case '.tif' 

137 tiftoload=filename; 

138 case '.csv' 

139 filetoload=filename; 

140 otherwise 

141 try 

142 filetoload=filename; 
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143 catch 

144 errordlg(lasterr,'File Type Error', 'modal') 
145 end 

146 end 

147 end 

148 end 

149 3 





150 % Read the current directory and sort the names 

151 § ------------------------------------------------------------ 
152 function load_listbox(dir_path, handles) 

153 cd (dir_path) 

154 dir_struct = dir(dir_path); 


155 [sorted_names,sorted_index] = sortrows({dir_struct.name}'); 
156 handles.file names = sorted_names; 
157 handles.is_ dir = [dir_struct.isdir]; 


158 handles.sorted_index = [sorted_index]; 
159 guidata (handles.newfirsttry, handles) 
160 set (handles.listbox1, 'String',handles.file_names,... 
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161 'value',1) 

162 set (handles.textl, 'String',pwd) 

163 

164 8 ——————— —S 

165 % --- Executes during object creation, after setting all properties. 
166 function listboxl_CreateFcen(hObject, eventdata, handles) 

167 % hObject handle to listboxl (see GCBO) 

168 % eventdata reserved - to be defined in a future version of MATLAB 
169 % handles empty - handles not created until after all CreateFcns called 
170 

171 $% Hint: listbox controls usually have a white background, change 
172 & "‘usewhitebg' to 0 to use default. See ISPC and COMPUTER. 


173 usewhitebg = 1; 
174 if usewhitebg 








7S set (hObject, 'BackgroundColor', 'white'); 

176 else 

a A i set (hObject, 'BackgroundColor', get (0, 'defaultUicontrolBackgroundColor')); 
178 end 

179 

180 3 _ 

181 % --- Executes during object creation, after setting all properties. 


182 function edit2_CreateFcn(hObject, eventdata, handles) 
183 % hObject handle to edit2 (see GCBO) 


184 % eventdata reserved - to be defined in a future version of MATLAB 

185 % handles empty - handles not created until after all CreateFcns called 
186 

187 % Hint: edit controls usually have a white background on Windows. 

188 3% See ISPC and COMPUTER. 
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189 
190 
191 
192 
193 
194 
195 


196 
197 
198 
199 
200 
201 
202 


203 


204 
205 
206 
207 
208 
209 
210 
211 
212 
213 
214 
215 


216 


aa 
218 
219 
220 


221 


222 
223 
224 
225 
226 
227 


228 
229 


if ispc 

set (hObject, 'BackgroundColor', 'white'); 
else 

set (hObject, 'BackgroundColor',get(0, 'defaultUicontrolBackgroundColor')); 
end 


op 
K 








% --- Executes on mouse press over axes background. 

function newname_ButtonDownFcn (hObject, eventdata, handles) 

% hObject handle to newname (see GCBO) 

$ eventdata reserved - to be defined in a future version of MATLAB 
$ handles structure with handles and user data (see GUIDATA) 


K 





% --- Executes on button press in pushbuttonl. (Vector Data Manipulator = vdatamanipuy” 
lator) 
function pushbuttonl_Callback (hObject, eventdata, handles) 
% hObject handle to pushbuttonl (see GCBO) 
$ eventdata reserved - to be defined in a future version of MATLAB 
$ handles structure with handles and user data (see GUIDATA) 
global filetoload aflag gflag foldername fitflag; 
if (aflag==[]) aflag=0; 
end 
if (gflag==[]) gflag=0; 
end 
if (fitflag==([]) fitflag=0; 
end 
Determine if this is a current session or a return to a previous’ 
session of data creation 
button = questdig('Are you in an active session of data analysis?','Session Type Selew 
ction'); 
switch button 
case 'Yes' 
nott='not'; 
[Valmanac, xnplot, ynplot,xinterp, yinterp, sample name, BeamEnergy,probe current] # 
=vdatamanipulator(filetoload, aflag, gflag, foldername) ; 
almanac_save(Valmanac, foldername, xnplot, ynplot, xinterp, yinterp, gflag,nott,filw 


etoload); 
Lf £itflag 
datatofit=[xinterp; yinterp]'; 
fbsLd(datatofit, foldername,sample_ name, BeamEnergy, probe current) ; 
end 
case 'No' 


prompt = {'Enter the folder name for data plot storage','Enter name of superpw 
osition plot to amend'}; 

dlg_ title = 'Amended Session Input'; 

num_lines= 1; 
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230 def Sty oe 

231 answer = inputdlg(prompt,dlg_title,num_lines, def); 

232 answerstr=char (answer) ; 

233 of=answer (1); 

234 oldfolder=char (of) 

235 superplotname=char (answerstr(2,:)); 

236 (Valmanac, xnplot, ynplot,xinterp, yinterp, sample name, BeamEnergy,probe current] @ 
=vdatamanipulator(filetoload, aflagl, gflagl, oldfolder); 

23) almanac_save(Valmanac, oldfolder, xnplot, ynplot, xinterp, yinterp, gflag, superplot ¢ 
name, filetoload) ; 

238 if fitflag 

239 datatofit=[xinterp;yinterp]'; 

240 fbsLd(datatofit, foldername, sample name, BeamEnergy, probe current); 

241 end 

242 otherwise 

243 disp('To find the superplot and folder to amend, look in SEM\. Folder data sev 
ts are named by a lxLetter designator denoting sample, followed by the beam energy. d 
The superplot is in this folder.') 

244 end 

245 

246 3% D4 

247 % --- Executes on button press in pushbutton2 (ImageDataManipulator) 

248 function pushbutton2_Callback(hObject, eventdata, handles) 

249 % hObject handle to pushbutton2 (see GCBO) 

250 $% eventdata reserved - to be defined in a future version of MATLAB 

251 $ handles structure with handles and user data (see GUIDATA) 

252 global tiftoload foldername aflag gflag fitflag; 

253 if (aflag==[]) aflag=0; 

254 end 

255 if (gflag==[]) gflag=0; 

256 end 

257 if (fitflag==[]) fitflag=0; 

258 end 

259 %Determine if this is a current session or a retué 
rn to a previous session of data creation 

260 button = questdlg('Are you in an active session of data analysis?','Session Type Selew 
ction', 'Cancel'); 

261 switch button 

262 case 'Yes' 

263 nott='not'; 

264 {almanacR, almanacC,xnplot, ynplot,xinterp, yinterp, sample name, BeamEnergy, prob ¢ 
e_current]=imagedatamanipulator(tiftoload, aflag, gflag, fitflag, foldername) ; 

265 almanac_save(almanacR, foldername, xnplot, ynplot, xinterp, yinterp,gflag,nott,tif @ 
toload) ; 

266 flagc=0; $Sets the superimposition flag to '0' so column data # 
will not be printed on superposition plot 

267 almanac_save(almanacC, foldername, xnplot, ynplot, xinterp, yinterp, flagc,nott,tif 


toload) ; 
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268 if fitflag 

269 datatofit=[xinterp;yinterp]'; 

270 fbsLd(datatofit, foldername, sample name, BeamEnergy, probe current); 

271 end 

272 case 'No!' 

273 prompt = {'Enter the folder name for data plot storage','Enter name of superpY# 
osition plot to amend'}; 

274 dig_title = 'Amended Session Input'; 

275 num_lines= 1; 

276 def Se fey nye 

277 answer = inputdlg(prompt,dlg_ title,num_lines, def) ; 

278 answerstr=char (answer) ; 

279 oldfolder=(answerstr(1,:)); 

280 superplotname=(answerstr(2,:)); 

281 nott='not'; %Ensures that column data will not be 

282 {almanacR, almanacC, xnplot, ynplot,xinterp, yinterp, sample name, BeamEnergy, prob @ 
e_current]=imagedatamanipulator(tiftoload, aflag, gflag,fitflag, oldfolder); 

283 almanac_save(almanacRk, foldername, xnplot, ynplot,xinterp, yinterp, gflag, superplo ¢ 
tname, tiftoload) ; 

284 flagc=0; $Sets the superimposition flag to '0' so column data will” 
not be printed on superposition plot 

285 almanac_save(almanacc, foldername, xnplot, ynplot,xinterp, yinterp, flagc,nott,tif # 
toload) ; 

286 if fitflag 

287 datatofit=[xinterp;yinterp]'; 

288 fbsLd(datatofit, foldername, sample name, BeamEnergy,probe current); 

289 end 

290 otherwise 

291 disp('To find the superplot and folder to amend, look in SEM\. Folder data sets arw# 
e named by a lxLetter designator denoting sample, followed by the beam energy. The sw 
uperplot is in this folder.') 

292 

293 end 

294 

295: §_ ame ee. 

296 % --- Executes on upcheck of checkboxl = This is for 'Keep Almanac for mulitple filesw# 
' 

297 function checkboxl_Callback(hObject, eventdata, handles) 

298 % hObject handle to checkboxl (see GCBO) 

299 % eventdata reserved - to be defined in a future version of MATLAB 

300 % handles structure with handles and user data (see GUIDATA) 

301 % Hint: get(hObject,'Value') returns toggle state of checkboxl 

302 global aflag; 

303 aflag = 1; 

304 

305 $2 ie ——— en 2 

306 % --- Executes on upcheck of checkbox2 for superimposition box. 
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348 {[legend_h, object_h,plot_h,text_strings] = legend(gca) ; 
349 [n,m]=size(text_strings); 

350 counter=m/2; 

35% legendname=text_ strings; 

352 if ishold 

353 yes=1; 

354 else hold on; 

355 end 

356 firsttime=0; 

357 end 


358 fid = fopen([{'C:\Documents and Settings\FMBradley\My Documents\Physics\Data Analysis\¥” 
', folder, '\Almanac.csv'],'a'); 

359 fprintf(fid,'$14.12f %d $d %6.3f %d %6.3f\n',Dalmanac) ; 

360 fclose(fid); 


361 if firsttime== print superimposed graphs of row vector data 
362 legendname={}; 

363 legendname=char {legendname) ; 

364 legendname=cellstr (legendname) ; 
365 counter=1; 

366 hold on; 

367 grid on; 

368 firsttime=0; 

369 end 

370 if flag == 

Sil switch counter 

372 case {1} 

oi plot (xnplot, ynplot,'og'); 
374 plot (xinterp, yinterp,'.g'); 
375 case {2} 

376 plot (xnplot, ynplot, 'ob'); 
377 plot (xinterp, yinterp,'.b'); 
378 case {3} 

379 plot (xnplot, ynplot,'oc'); 
380 plot (xinterp, yinterp,'.c'); 
381 case {4} 

382 plot (xnplot, ynplot, 'ok'); 
383 plot (xinterp, yinterp,'.k"'); 
3B4 case {5} 

385 plot (xnplot, ynplot,'om'); 
386 plot (xinterp, yinterp,'.m'); 
387 case {6} 

388 plot (xnplot, ynplot, 'oy'); 
389 plot (xinterp, yinterp,'.y'); 
390 case {7} 

391 plot (xnplot, ynplot,'or'); 
392 plot (xinterp, yinterp,'.r'); 
393 case {8} 

394 plot (xnplot, ynplot,'sg'); 
395 plot (xinterp, yinterp,'*g'); 
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41 end 

42 %Test Row data for least variance of column numbers for selected column data sample 

43  I0T=10'; 

44 (Y2,13]=max(I0T); $Two vectors: Y2=max value of each Row, I3=Col# of each row's #@ 
Max) 

45  [KR,LR]=size (13); 

46 

47 MinVarianceR=500; 

48 a=25; 

49 for m3=2: (LR-a-1) $Loop through all row numbers 

50 for n=l:a Begin creation of test vector. Populate with 25 elements of I¥ 
2 (Column # of max from each row) 

51 TestR(n)=13((m3+n)); 

52 end 

53 VarianceR=std(TestR) ; $Calculate variance of maximum values' column #s to #¥ 
select desired data sample. 

54 if VarianceR < MinVarianceR $Update minimum variance 

55 MinVarianceR=VarianceR; 

56 seqgstartR=m3; $save indice start point in vector I2 of desired dag@¢ 
ta sample (Column #s of desired data sample) 

2 seqstopR=m3+a; 

58 seqlengthR=a; define sequence length to allow for seqlength to bew 
used as indice addition term for creation of new vector 

59 end 

60 end 

61 

62 %Compare MaxValue vectors and choose column and row with largest same maximum 

63 PeakPixelValue=0; 

64 for stepR=1:25; 

65 for stepC=1:25; 

66 if Y(stepC+seqstart) ==Y2 (stepR+seqstartR) 

67 Peak=Y(stepC+seqstart) ; 

68 if Peak>PeakPixelValue 

69 PeakPixelValue=Peak; 

70 MaxPixelColNum=seqstart+stepC; 

vee MaxPixelRowNum=seqstartR+stepR; 

72 end 

73 end 

74 end 

75 end 

76 

77 extract row and column data 

78 RowData=10 (MaxPixelRowNum,:); 

79 ColData=I0(:,MaxPixelColNum) ; 

80 

81 %Create Noise vector from data outside of spike 

82 %Calculate variance of noise and through out sample data vectors 

83 for z=1:50 

84 NoiseData (z)=RowData(z); 
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85 end 
86 
87 RawVarOfNoise=std(NoiseData) ; 
88 RawMeanOfNoise=mean(NoiseData); 
89 for z1=1:50 
90 if NoiseData(zl) >= RawMeanOfNoise+RawVarOfNoise 
91 NoiseData(z1)=RawMeanOfNoise; 
92 end 
93 end 
94 MeanOfNoise=mean (NoiseData); $Calculate average of noise within variance 
95 vVarOfNoise=std(NoiseData) ; 
96 $Normalize row vector of sample data 
97 NormRowData=(RowData-MeanOfNoise) / (PeakPixelValue-MeanOfNoise) ; 
98 
99 Normalize ColData vector 
100 for z=1:50 
101 NoiseDataCol (z)=ColData(z); 
102 end 
103 
104 RawVarOfNoiseCol=std(NoiseDataCol); 
105 RawMeanOfNoiseCol=mean (NoiseDataCol) ; 
106 for z1=1:50 
107 if NoiseDataCol(z1) >= RawMeanOfNoiseCol+ (2*RawVarOfNoiseCol) 
108 NoiseDataCol (z1)=RawMeanOfNoiseCol; 
109 end 
110 end 
a Ms 
112 %Calculate average of noise within variance 
113 MeanOfNoiseCol=mean(NoiseDataCol); 
114 varOfNoiseCol=std(NoiseDataCol) ; 
115 WNormColData=(ColData-MeanOfNoiseCol) / (PeakPixelValue-MeanOfNoiseCol) ; 
116 
117 ‘%add spline interpolation 
118 xr= (MaxPixelColNum-49) : (MaxPixelColNum+50) ; 
119 xc= (MaxPixelRowNum-49) : (MaxPixelRowNum+50) ; 
120 for x1=1:100 
127 yr(xl) = NormRowData (MaxPixelColNum-50+x1) ; 
ny yc (x1) =NormColData (MaxPixelRowNum-50+x1) ; 
123 end 
124 csc = spline(xc,[0 ye 0]); 
125 csr=spline(xr,[0 yr 0]); 
126 xxr=linspace ( (MaxPixelColNum-50), (MaxPixelColNum+50),1000); 
127 xxc = linspace( (MaxPixelRowNum-50), (MaxPixelRowNum+50),1000); 
128 
129 
130 $%Extract Halfmaxfullwidth from the normalized Row data 
131 maxdif=1;maxdify=1; 
132 for xxx=1:500 
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133 dif=abs (0.2-ppval (csr,xxr(xxx))); 
134 if dif<maxdif 
135 maxdif=dif; 
136 lhshalfmax=xxr (xxx); 
137 end 
138 end 
139 for xxy=500:1000 
140 dify=abs (0.2-ppval (csr, xxr(xxy))); 
141 if dify<maxdify 
142 maxdify=dify; 
143 rhshalfmax=xxr (xxy); 
144 end 
145 end 
146 hold; 
147 = Rhalfmaxfullwidth=rhshalfmax-lhshalfmax 
148 
149 %Extract the halfmaxfullwidth from the normalized column data 
150 maxdif=.2;maxdify=.2; 
151 for xxx=1:500 
152 dif=abs (0.2-ppval(csc,xxc (xxx) )); 
153 if dif<maxdif 
154 maxdif=dif; 
155 lhshalfmax=xxc (xxx) ; 
156 end 
157 end 
158 for xxy=500:1000 
159 dify=abs (0.2-ppval (csc, xxc (xxy))); 
160 if dify<maxdify 
161 maxdify=dify; 
162 rhshalfmax=xxc (xxy) ; 
163 end 
164 end 
165 hold; 
166 Chalfmaxfullwidth=rhshalfmax-lhshalfmax; 
167 
168 %Assign output variables to almanac vector 
169 prompt = {'Unique sample name','Enter probe current in "3e-8" notation:','Enter Beam #@ 
Voltage in format "30" kV:', 'Enter Exposure Time in "10.005" (sec) format:'}; 
170 dilg_title = 'Almanac File Definition Input'; 
171 num_lines= 1; 
TZ deft = {"'S72",'6e-10', °25°, °1.0*}7 
173 answer = inputdlg(prompt,dlg_title,num_lines, def); 
174 answerstr=char (answer) ; 
175 sample name=(answerstr(1,:)); 
176 probe_current=(answerstr(2,:)); 
177 BeamEnergystring=(answerstr(3,:)); 
178 BeamEnergy=str2double(answerstr(3,:)); 
179 ExpTime=str2double(answerstr(4,:)); 
180 almanacR=[sample_ name str2double (probe current) 0 BeamEnergy Rhalfmaxfullwidth PeakPiw 
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181 


182 
183 
184 
185 
186 
187 
188 
189 


190 
191 
192 
193 
194 
195 
196 
197 
198 
199 
200 
201 
202 
203 
204 
205 
206 
207 
208 
209 
210 
214 
212 
213 
214 
215 
216 
217 
218 


219 
220 
221 
222 


xelValue ExpTime]; 
almanacC=[sample_ name str2double (probe current) 1 BeamEnergy Chalfmaxfullwidth PeakPil¥ 
xelValue ExpTime]; 


% Page 1 graphs. Overview of image with selected max values in colmax/row 
$ and rowmax/col 


figure (1) 
title(['Overview page for: ',folder,'_',tiftoload]); 
subplot (3,1,1), imagesc(abs(I0)); $Shows reproduction of tiff f¥” 


ile from MicroCCD 

hold; 

plot (MaxPixelColNum, MaxPixelRowNum, 'xk'); 
title(['MATLAB reproduction of TIFF (',tiftoload,')']); 
hold off; 

$ Plots maximum per Row vs Column number where that maximum falls in the row. 
subplot (3,1,2), plot(max(I0)); 

hold; 

title(['Max Pixel value/Row vs Column Number of MaxPixVal (',tiftoload,')']); 
plot (MaxPixelColNum, (MeanOfNoise-VarOfNoise),'*y'); 

axis ([0,L, (MeanOfNoise-VarOfNoise), (max(Y¥)+10)]); 

grid on; 

hold; 


%$ Graph Columns' max pixel value vs row Number of MaxpixVal 

subplot (3,1,3), plot (max(IOT)); 

hold; 

title(['Max Pixel value/Column vs Row Number of MaxPixVal (',tiftoload,')']); 
plot (MaxPixelRowNum, (MeanOfNoiseCol-VarOfNoiseCol),'*y'); 

axis([0,LR, (MeanOfNoiseCol-VarOfNoiseCol), (max (¥2)+10)]); 

grid on; 

hold; 


$Save figure 1 to file 

D=now; 

ddate=day (D); 

{n,mmonth]=month(D); 

dateserial=['_',num2str(ddate) ,mmonth, num2str(fileserial) ]; 

filetosave=['C:\Documents and Settings\FMBradley\My Documents\Physics\Data Analysis\' ¥¢ 
,folder,'\OV',sample_name,'_',num2str(BeamEnergy),'_',probe_ current, dateserial,'.fig'¥ 
l; 


saveas (1, filetosave); 


$Write NormRowData to file 

Vdatafile=['C:\Documents and Settings\FMBradley\My Documents\Physics\Data Analysis\', 4 
folder, '\VRdata_',sample_name,'_',num2str(BeamEnergy),'_',probe_current,dateserial,'. ¢ 
esv'" 1; 
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223 fid = fopen(Vdatafile,'w'); 

224 fprintf(fid, '$11.9f£\n',NormRowData) ; 

225 fclose(fid); 

226 

227 Write NormColData to file 

228 Vdatafile=['C:\Documents and Settings\FMBradley\My Documents\Physics\Data Analysis\', # 
folder, '\VCdata_',sample_name,'_',num2str(BeamEnergy),'_',probe_current,dateserial,'.” 





esv']; 
229 fid = fopen(Vdatafile,'w'); 
230 fprintf(fid, '$11.9f\n',NormColData) ; 
231 f£fclose(fid); 
232 
233 %Plot Raw row and column data sets 
234 figure(2); 
235 title(['Peak plots for: ',folder,'_',tiftoload]); 
236 %Raw Row and Column data 
237 subplot(2,2,1), plot (RowData,'o'); 
238 title(['Raw Data for Peak Pixel Row (',tiftoload,')']); 
239 axis([(MaxPixelColNum-50), (MaxPixelColNum+50),0, (PeakPixelValue+10)]); 
240 grid on; 
241  subplot(2,2,2), plot(ColData,'o'); 
242 title(['Col Data for Peak Pixel Column (',tiftoload,')']); 
243  axis([(MaxPixelRowNum-50), (MaxPixelRowNumt+50),0, (PeakPixelValue+20)]); 
244 grid on; 
245 
246 %Print Normalized Row Data 
247 R=0;xRlimit=0;xRstep=0; 
248 [R,xRlimit]=size(NormRowData) ; 
249 for xRstep=1:xRlimit; 


250 XRnorm(xRstep) =xRstep-MaxPixelColNum; 
251 end 

252 subplot (2,2,3), plot(XRnorm, NormRowData, 'o'); 
253 hold; 


254 plot ((xxr-MaxPixelColNum),ppval(csr,xxr),'.'); 
255 title(['Normalized Row Data for Peak Pixel Row (',tiftoload,')']); 
256 axis([(-50),(50),-.2,1.1]); 

257 grid on; 

258 text(10,0.5, ['FWHM=',num2str (Rhalfmaxfullwidth) ]); 
259 

260 

261 

262 

263 %Print Column Data 

264 [xlimit,C]=size (NormColData) ; 

265 for xstep=1:xlimit; 


266 XCnorm(xstep) =xstep-MaxPixelRowNum; 

267 end 

268 subplot(2,2,4), plot(XCnorm, NormColData,'o'); 
269 hold; 
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270 = plot((xxc-MaxPixelRowNum) ,ppval(csc,xxc),'.'); 

271 title(['Normalized Col Data for Peak Pixel Column (',tiftoload,')']); 

272 axis({(-50),(50),-.2,1.1]); 

273 grid on; 

274 text(10,0.5, ['FWHM=',num2str (Chalfmaxfullwidth) ]); 

27D 

276 %Save figure 2 to file 

277 filetosave2=['C:\Documents and Settings\FMBradley\My Documents\Physics\Data Analysis\¥ 
', folder, '\RC_',sample_name,'_',num2str(BeamEnergy),'_',probe_ current, dateserial, '.¢ 
fig']; 

278 saveas(2, filetosave2); 

279 

280 $Test for desire for superimposed graph page and set export variables for 

281 %consolidation in newfirsttry 

282 if (or(gflag,fitflag) ) 

283 xnplot=XRnorm; 

284 ynplot=NormRowData; 

285 xinterp=xxr-MaxPixelColNum; 

286 yinterp=ppval (csr,xxr); 

287 else 

288 xnplot=0; 

289 ynplot=0; 

290 xinterp=0; 

291 yinterp=0; 

292 end 
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C:\Documents and Settings\FMBradley\My Documents\Phy...\fbsLd.m Page 1 
November 16, 2005 5:43:45 PM 
1 . Created by FM Bradley to perform a least squares fit for data to the 
2 $2-D model for diffusion of minority carriers in a luminscent 
3  %semiconductor. 
S 
5 Revision Log: 13 September - Time remaining window added by FMBradley 
6 
7 $Takes a 2-D vector data in (xstepX2) size, normalizes, cuts to 1% of max, 
8 %performs a least squares fit to a 2-D model of carrier distro given a 
9 %guassian distributed generation area. Output is the fitted model vector, 
10 %the difference between model and experimental data, plots of the parameter 


wh tspace, both curves, the selected of optimal Mu, and n=2sigma radius of 
1Z generation area that fits the data. 


14 function [FittedCurve, residue, Paramspace, Lopt,nopt]=fbsLd(datain, folder,sample name,BMl 
eamEnergy,probe current) 











15 

16 $%CreateStruct.WindowStyle='replace'; 

17 $ CreateStruct.Interpreter='tex'; 

18 "Calculating: 3:6.<.4 ', 'Progress Monitor',CreateStruct) ; 

19 nmin=le-4 $Significant radius of generation volume distribution(cm) 95% of currew# 
nt ring?? 

20 v=0; $Applied bias (Volts) 

21 tau=4.32e-9; $input variable. Lifetime of carrier in question (s) 

22 icd=1.0; sinter-contact distance (mm) for calculation of electric field assoc¥ 
iated with V 

23 E=V/icd*10; Electric field (V/cm) 

24  step=0.04e-4; $step size in (cm)designed to fit the incré@ 
ement size of our experimental data w/interpolation 

25 [R3,13]=size(datain) ; $Id dimension of datain array 

26 datain(:,1)=0.4e-4.*datain(:,1); $Use this line if the data is coming from w# 
actual experimentally gathered and spline interpreted data from the MATLAB GUI. 

27 ‘Necessary because spline data is already 

in microns and we must convert all numbers here to cm "4 

28 Lmin=SlopeL2 (datain)-.3e-4; sConducts in slope analysis method for t#@¢ 
he diffusion length assessment 

29 Lstep=.le-4; 

30 stop=0; 

31 hile (stop<3) 

32 switch stop 

33 case 0 

34 Lstep=.2e-4; 

35 nstep=.2e-4; 

36 Lmax=Lmin+.4e-4; 

3 nmax=nmin+.4e-4; 

38 case 1 

39 Lstep=.1le-4; 

40 nstep=.le-4; 

41 Lmax=Lmin+.2e-4; 
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42 nmax=nmin+.2e-4; 
43 case 2 
44 Lstep=.le-4; 
45 nstep=.le-4; 
46 nmin=nmin-.5e-4; 
47 Lmin=Lmin-.5e-4; 
48 Lmax=Lmin+le-4; 
49 nmax=nmin+le-4; 
50 end 
5. 
52 a=0; $step variable for indexing x 
53 b=0; $step variable for indexing y 
54 MinSum=10000; 
55 istep=0; 
56 Paramspace=0; 
S57 for L=Lmin:Lstep:Lmax 
58 t0=cputime; 
59 n=nmin; $1-d spot dimension of a squag¢ 
re generation area (cm) 
60 istep=istep+1; $parameter space indexing var #@ 
iable for mu 
61 jstep=0; $param space index variable fv 
orn 
62 S=(E*L*2) /.025; Field effect (cm) 
63 mtproduct=L*2/.025; $mu tau product in (cm*2/V) 
64 for n=nmin:nstep:nmax 
65 G=2/n; Gaussian 2sigma radius f#¢ 
or 95% generation distro (unitless) 
66 jstep=jsteptl; 
67 xmin=datain(1,1); convert datain from pixels to cm dew 
pendent upon resolution size of SEM/Optical scope .4microns/pixel 
68 xmax=max (datain(:,1)); 
69 x=xmin; "4 
% (cm) 
70 ymin=0; "4 
% (cm) 
pe: ymax=0; 4 
$ (cm) 
BA y=ymin; 
73 b=0; 
74 for yl=ymin:step: ymax 
75 a=0; 
76 b=b+1; 
te for xl=xmin:step:xmax 
78 Integrand=0; 
79 a=a+l; 
80 Integrand=dblquad(@Intgrnd, (xl-n), (xl+n), (yl-n), (yl+n),step,[],xl,yl, @ 
S,L,G); 
81 if (isnan(Integrand)) Integrand=OldIntegrand+.01*OldIntegrand; 
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82 end 
83 OldiIntegrand=Integrand; 
84 Int (a,b,1)=x1; 
85 Int (a,b, 2)=yl; 
86 Int(a,b,3) = (G/(2*pi*2*L*2) ) *Integrand; 
87 end send x loop 
88 end tend y loop 
89 %Least square computation 
90 (M2,12]=max(Int(:,1,3)); 
91 NormInt=[Int(:,1,1) (Int(:,1,3)./M2)]; 
92 $LogNormInt=log(Normint (:,2)); 
93 $LogDataIn=log (normdatain(:,2)); 
94 LS=0; 
95 
96 for loopc=1:R3; $1: (((dimofdata-1) /2)-105); Z 
$start loop at max +13 and count to end testing for least squares diff 
97 $LS (loopc) =(LogNormInt (I12+100+loopc) -LogDataIn(13+100+loopc) ) *2 
98 $; Used to calculate log difference squares between two lines 
99 Sof data 
100 if isnan(NormInt (loopc, 2) ) 
1O1 NormInt (loopc, 2)=1; 
102 end 
103 if isnan(datain(loopc, 2) ) 
104 datain(loopc, 2)=1; 
105 end 
106 if datain(loopc,2)<.7 Eliminates t¢ 
he top 30% of the curve that does not fit our assumptions 
107 LS (loopc) = (Normint (loopc, 2)-datain(loopc, 2) )*2; 
108 end ends non mod#@ 
el fit data 
109 end tends Least s¥ 
quares sum loop 
110 LoopSum=sum (LS) ; 
111 Paramspace (istep, jstep, 1)=(istep-1) *Lstep+Lmin; 
112 Paramspace (istep, jstep, 2)=(jstep-1) *nstep+nmin; 
113 Paramspace (istep, jstep, 3) =LoopSum; 
114 if LoopSum<MinSum 
115 MinSum=LoopSum; 
116 residue=MinSum; 
117 ModelFit=NormInt; 
118 $Muopt=mu; 
119 nopt=n; 
120 Lopt=L; 
t22 imin=istep; 
122 jmin=jstep; 
123 end 
124 figure(7); $Plot the search space each round of calcw 


ulations complete 
125 nplot=n*le4; 
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126 Lplot=L*1e4; 

127 if (stop==1)plot (Lplot,nplot,'xr'); 

128 else plot(Lplot,nplot, 'xk'); 

129 end 

130 if not(ishold) hold; 

132 end 

132 end tend of n for loop 

133 % CreateStruct.WindowStyle='replace'; 

134 % CreateStruct.Interpreter='tex'; 

135 % h=msgbox(['Total time remaining= ',num2str(((cputime-t0) /60)*(Lmax-L)/Lstep),' m#@ 
inutes. Step ',int2str(istep*jstep),' of ',num2str((((Lmax-Lmin) /Lstep) +1) *(((nmax-nm#¢ 
in) /nstep)+1)),' Total Steps'],'Progress Monitor',CreateStruct) ; 

136 end tend of L for loop 

137 

138 %Decision tree for creating direction of propogation of parameter space 

139 %build up. 

140 switch imin 

141 case 1 

142 Lmin=Lmin-Lstep; 

143 if jmin== 

144 nmin=nmin-nstep; 

145 elseif jmin== 

146 nmin=nmin; 

147 else nmin=nmin+nstep; 

148 end 

149 case 2 

150 Lmin=Lmin; 

151 if jmin== 

152 nmin=nmin-nstep; 

153 elseif jmin==3 

154 nmin=nmin+nstep; 

155 else 

156 stop=stopt1; 

L57 plot ( (Lmin*1le4), (nmin*le4), 'xb"'); 

158 end 

159 case 3 

160 Lmin=Lmin+Lstep; 

161 if jmin== 

162 nmin=nmin-nstep; 

163 elseif jmin==2 

164 nmin=nmin; 

165 else nmin=nmintnstep; 

166 end 

167 otherwise 

168 break; 

169 end tend for the switch loop 

170 end $end for the while loop 

Li1 title('Search Space'); Label Figure 7 

72 xlabel('Ld (\mum)'); 
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WIS ylabel('2\sigma (\mum)'); 

174 hold; Hold off for figure 7 

175 FittedCurve=ModelFit; 

176 %Write Fitted Curve data to file 

177 Vdatafile=['C:\Documents and Settings\FMBradley\My Documents\Physics\Data Analysis\', # 
folder, '\ModelFit',sample_ name, 'Ld_',Lopt,'_',num2str(BeamEnergy),'kV_',probe_current # 
5 cee 

178 fid = fopen(Vdatafile,'w'); 

179 fprintf(fid,'%11.9f\n',ModelFit) ; 

180 £fclose(fid); 

181 Residue=MinSum 

182 nopt 

183 Lopt 

184 

185 

186 %plot the resulting fitted curves 

187 figure(8); 

188 plot (ModelFit(:,1),ModelFit(:,2),'.b',datain(:,1),datain(:,2),'-r'); 

189 hold on; 

190 Xmax=xmax; 

191 Xmin=xmin; 

192 ymin=0; 

193 yYmax=1.0; 

194 Axis([Xmin Xmax Ymin Ymax]); 

195 title('Normalized Experimental vs. Model Fit"); 

196 xlabel('Radial distance from beam center (cm)'); 

197 ylabel('Normalized Intensity'); 

198 = $"\mum: \mu=', num2str (Muopt), ':\tau=',num2str(tau*le9), 

199 legendnamel=['Model Fit: (n=',num2str(nopt*le4),'\mum: Ld=',num2str(Lopt*le4),'\mum) Rw 
esidue=',num2str (residue) ]; 

200 legendname2=['Data: ',sample name,':PC=',num2str(probe current), 'Amps']; 

201 legend(legendnamel, legendname2) ; 

202 hold off; 

203 

204 

205 %Plot the parameter space of the resulting 1MicronxXlMicron parameter range for both 2¥¢ 
sigma and Ld 

206 figure(9); 

207 contour (Paramspace(:,:,2),Paramspace(:,:,1),Paramspace(:,:,3),100); 

208 hold on; 

209 plot(nopt,Lopt, '*y') 

210 Xmax=nmax; 

211 Xmin=nmin; 

212 Ymin=Lmin; 

213  Ymax=Lmax; 

214 Axis([Xmin Xmax Ymin Ymax]); 

215 title([{'Parameter Space: ',sample_name,':',num2str(BeamEnergy),'kV:',num2str(probe cu 
rrent),'A']); 

216 ylabel('Ld step multiples [cm*2/Vs]'); 
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APPENDIX D. VECTOR DATA EXTRACTION ROUTINE 
“VDATAMANIPULATOR.M” 





N 





3 to 4 
of multiple calls to 
: s are beingwv 
to file, 





and increments fileserial 





6 to help differeniate 


es 


rh 
He 
.Q 
ty 


name:','s'); 















IMM NMNMNDND Fe 








MeanOfNoise=me Calcula 
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47 
48 
49 
50 
S51 
52 
53 
54 
55 
56 
7 
58 
99 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 


70 
71 
72 
73 
74 
25 
76 
va 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 


89 
90 
91 
92 


csr=spline(xr,[0 yr 0]); 


Page 2 
5:46:14 PM 


xxr=linspace ( (MaxPixelColNum-50), (MaxPixelColNum+50),1000); 


$Extract Halfmaxfullwidth from the normalized Row data 


maxdif=1;maxdify=1; 
for xxx=1:500 


dif=abs (0.2-ppval(csr,xxr(xxx))); 


if dif<maxdif 
maxdif=dif; 
lhshalfmax=xxr (xxx) ; 
end 
end 
for xxy=500:1000 


dify=abs (0.2-ppval (csr, xxr(xxy)));7 


if dify<maxdify 
maxdify=dify; 
rhshalfmax=xxr(xxy); 
end 
end 


Rhalfmaxfullwidth=rhshalfmax-lhshalfmax; 


$Assign output variables to almanac vector 
prompt = {'Unique sample name','Enter probe current 


in "3e-8" notation:','Enter Beam #¢ 


Voltage in format "30" kv:','Enter Exposure Time in "10.005" (sec) format:','Enter R # 


or C for Row or Column Data'}; 


dlg_title = 'Almanac File Definition Input'; 


num_lines= 1; 


def = {"S72" ,.16e=10',.'255, '1.0' "RY he 
answer = inputdlg(prompt,dlg_title,num_lines,def) ; 


answerstr=char (answer) ; 

sample name=(answerstr(1,:)); 
probe current=(answerstr(2,:)); 
BeamEnergystring=(answerstr(3,:)); 


BeamEnergy=str2double(answerstr(3,:));7 


ExpTime=str2double(answerstr(4,:));7 
switch answerstr(5) 
case 'R' 
vector type='0'; 
case 'C' 
vector type='1'; 
otherwise 


dialg('Your test did not work'); 


end 


almanacV=[str2double (probe current) vector_type BeamEnergy Rhalfmaxfullwidth PeakPixe @ 


lValue ExpTime]; 


$Analysis output loop 
$Plot Raw row and column data sets 
figure (1); 
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93 $%Plot overview graph of row 
94  subplot(1,2,1), plot(RowData,'o'); 
95 title(['Raw Data for Peak Pixel Vector (',foldername, filetoload,')']); 
96 axis([(MaxPixelColNum-50), (MaxPixelColNum+50),0, (PeakPixelValue+10)]); 
97 grid on; 
98 
99 %Print Row Data 

100 [xRlimit, R]=size (NormRowData) ; 

101 for xRstep=1:xRlimit; 

102 XRnorm(xRstep) =xRstep-MaxPixelColNum; 

103 end 

104 _~ subplot(1,2,2), plot (XRnorm, NormRowData, '0'); 

105 hold; 

106 plot( (xxr-MaxPixelColNum) ,ppval(csr,xxr),'.'); 

107 title(['Normalized Data for Peak Pixel Vector (',filetoload,')']); 

108 axis([{(-50), (50) ,—.2,1.1]}; 

109 grid on; 

110 text(10,0.5, ['FWHM=',num2str (Rhalfmaxfullwidth) ]); 

111 = hold; 

112 

113 %Save figure 1 to file 

114 D=now; 

115 ddate=day(D); 

116 = [n,mmonth]=month (D) ; 

117 dateserial=[' ',num2str(ddate) ,mmonth, num2str(fileserial)]; 

118 filetosave=['C:\Documents and Settings\FMBradley\My Documents\Physics\Data Analysis\' ¢ 

,foldername, '\', answerstr(4),num2str(BeamEnergy),'_',probe_current,dateserial,'.fig'] ¢ 

119  saveas(1,filetosave) ; 

120 figure(2); 

E21 

122 %Set export variables for 

123 ‘%consolidation in newfirsttry 

124 xnplot=XRnorm; 

LAS ynplot=NormRowData; 

126 xinterp=xxr-MaxPixelColNum; 

27 yinterp=ppval (csr,xxr); 
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APPENDIX E. SLOPE ANALYSIS ALGORITHM (SLOPEL2 .M) 
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43, Lr=-1/pr(1) 

44 pr=pr (2) 

45 fitdataR=[xfitr', yfitr']; 

46 

A] %----- Plot ln(Intensity) vs X pos for all data left of spot------ 

48 figure(4); 

49 plot (datavector(:,1),log(datavector(:,2))); 

50 $%axis([(datavector(1,1)), (0), (log(min(datavector(:,2)))),(0)1); 

51 title('LN(Intensity) of Left Side Distribution'); 

52 xlabel('Radial distance from beam center (cm)'); 

53  ylabel('LN(Normalized Intensity) '); 

54 

55 $----- Calculate error reports and plot line slopes------ 

56 figure(5); 

57 Lplus=L+(.4e-4) / (max (yfit)-yfit(1)); $The Plus factor should be .4e-4 for 
real data and .04e-4 for tfbs control data 

58 Lminus=L-(.4e-4) /(max(yfit)-yfit(1)); $The minus factor should be .4e-4 forg 
real data and .04e-4 for tfbs control data 

59 plot(xfit,yfit,'ob'); 

60 hold; 

61 plot(xfit,polyval(p,xfit),'-r'); 

62 yplus=(1/Lplus) *xfit+b; 

63 yminus=(1/Lminus) *xfit+b; 

64 plot(xfit,yplus,'-g'); 

65 plot(xfit, yminus, '-k'); 

66 % Right side 

67 Lplusr=Lr+(.4e-4) / (max (yfitr)-min(yfitr) ); The Plus factor should be .4e-¢ 
4 for real data and .04e-4 for tfbs control data 

68 Lminusr=Lr-(.4e-4)/(max(yfitr)-min(yfitr)); $The minus factor should be .4e¥ 
-4 for real data and .04e-4 for tfbs control data 

69 plot(xfitr,yfitr,'ob'); 

70 = plot(xfitr,polyval(pr,xfitr),'-r'); 

71 =~ yplusr=(-1/Lplusr) *xfitr+br; 

72 = yminusr=(-1/Lminusr) *xfitr+br; 

73° =plot(xfitr,yplusr,'-g'); 

74° plot(xfitr, yminusr, '-k'); 

i fs) 

76 legendnamel=['Data Points']; 

77° legendname2=['L actual Ld=',num2str(L*le4),'\mum']; 

78 legendname3=['L Ld(+)=',num2str(Lplus*le4),'\mum']; 

79 legendname4=['L Ld(-)=',num2str (Lminus*1le4),'\mum']; 

80 legendname5=['Data Points']; 

81 legendname6=[('R actual Ld=',num2str(Lr*le4),'\mum']; 

82 legendname7=['R Ld(+)=',num2str (Lplusr*1le4),'\mum']; 

83 legendname8=['R Ld(-)=',num2str(Lminusr*1le4),'\mum']; 

84 

85 

86 legend(legendnamel, legendname2, legendname3, legendname4, legendname5, legendname6é, legend 


name7, legendname8) ; 


86 





el('Radial distance from beam center (cm)'); 


el('LN(Normalized Intensity) '); 
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APPENDIX F. NUMERICAL INTEGRATION SOLUTION FOR 
MINORITY CARRIER DISTRIBUTION (INTEGRAND .M) 
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